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The Role of Convergence in Stereoscopic Vision 
By W. DY WRIGHT 
Technical Optics Section, Imperial College, London 
Read before the Optical Group, 25th January 1950; MS. received 12th October 1950 


ABSTRACT. An apparatus for testing stereoscopic acuity has been constructed in which 
the distance of a reference aperture P is compared with that of a test aperture Q selected 
from one of a tapering series of aperture plates. The apertures are seen as illuminated discs 
in a black field and all monocular aids to depth perception are eliminated. The observer has 
to state whether Q is in front or behind P, and the acuity is derived by plotting the percentage 
estimates against the distance of Q@. Observations have been made for condition A, in which 
the observer is allowed to look first at P and then at Q, and for condition B, in which the 
observer fixates on P and views Q extra-foveally. From a comparison of the acuity for 
the two conditions for various angular separations 0 between P and Q, deductions can be 
made about the relative contributions to stereoscopic vision of convergence and retinal 
disparity. The results suggest that convergence makes a significant contribution and 
has become the predominant factor at 9=20°. For 0=14°, the acuity for condition A is 
due entirely to convergence, since the blind spot interferes with the retinal disparity. Tests 
with flash exposures suggest that small eye movements do not play any important part in 
depth perception. The presence of depth perception under conditions of gross diplopia was 
confirmed, indicating that fusion is not.an essential item for three-dimensional vision. 


Sf SIND ROD UCrLON 

ARIOUS theories of depth perception have been proposed from time to 

time and, in common with other aspects of vision, their popularity has 

waxed and waned. ‘The theories have not necessarily been mutually 
exclusive, since it is apparent that a number of factors can contribute to three- 
dimensional vision, yet sometimes one, sometimes another, has been regarded as 
pre-eminent. Dispute has chiefly centred on the relative importance of 
convergence and retinal disparity. Some authorities (for instance, Javal 1866) 
have held that convergence was the primary factor, a near object being recognized 
as such through the influence of the extra-ocular muscles as they were tensed to 
converge the eyes on to the object; for distant objects, the visual axes are more 
nearly parallel and the muscles relaxed. Briicke (1841) considered that the 
perception of depth arose from foveal viewing and fusion of successive points 
round the contours of an object. The view generally held at the present time is 
that retinal disparity makes the main contribution. Since the right and left eyes 
view a scene from slightly different positions, the images on the two retinae have 
slightly different geometrical configurations, that is, there will be retinal disparity. 
The mental fusion of these dissimilar images is then assumed to evoke the sense of 
depth (see, for example, Duke-Elder 1932, p. 1064). 
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The evidence which is generally used to clinch the disparity argument is 
derived from tests on stereoscopic acuity under flash illuminations (for example, 
Langlands 1926, 1929). Doubts about the interpretation of Langlands’ results, 
however, have been expressed (Wright 1948), and certainly in the case of extra- 
foveal vision, Langlands found no stereoscopic acuity under flash illuminations, 
Yet an excellent sense of depth can be appreciated between objects widely separated 
across the field of view, and it seems reasonable to suppose that this must largely 
result from the successive foveal convergence on the individual objects, since the 
poor resolving power of the peripheral retina makes it unlikely that retinal disparity 
in the outer parts of the field can make a very effective contribution to depth 
perception. 

When stereoscopic vision was discussed in 1948 at the Optical Convention 
organized by the Spectacle Makers Company, the need for fresh experiments to 
provide a more decisive test of the relative contributions of convergence and 
disparity was apparent. The matter is not merely of academic interest but has a 
practical bearing in connection with certain aspects of orthoptic training, and with 
problems of space perception under abnormal conditions, for example when the 
observer or the target is moving at high velocity, when the illumination level is low, 
etc. ‘The Spectacle Makers Convention has thus provided the stimulus for the 
present work. 

§2. APPARATUS 

In ordinary seeing, there are a number of subsidiary aids to the judgment of 
distance, such as size, perspective, parallax, etc., in addition to those derived 
from binocular vision. For the purpose of the present experiments, these aids 
had to be eliminated, and it was eventually arranged that the two objects whose 
distance was to be compared should consist of two illuminated, circular apertures 
seen against a black surround in a dark room. No structure of the surround 
could be resolved by the observer, hence he could gain no assistance in his judgment 
from the distance of the apertures relative to the background. This condition is 
difficult to satisfy except in tests with bright objects in a dark room, but this 
solution does impose some limitation on the state of adaptation under which the 
test can be applied. 


Aperture Plates; : 
a er 


wie Si ia} ae : : . Ag 
L IQ ae. 
2 | 6 ae ae Observer 
P ee ; 


Figure 1. Experimental arrangement for testing stereoscopic acuity. 


The reference object was provided by any one of four apertures P, Figure 1, 
located on an arc at a distance of 6 ft. from the observer. The test aperture Q was 
provided by any one of a battery of aperture plates mounted along a spindle so 
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that they could be swung into position in front of an illuminated diffusing 
screen L,. The apertures were spaced 1 in. apart over a distance of 14 in., except 
that intermediate plates at $ in. spacing were included in the central 5 in. The 
apertures formed a tapering series such that they all subtended the same angle 
(approximately 25 minutes of arc) when viewed by an observer 6 ft. from the 
mid-aperture. ‘The illuminated screen at the far end of the apertures provided a 
structureless area of light, against which Q appeared as a uniform disc of light 
defined only by its bounding edge. 

The observer’s head was located and fixed by biting on a dental impression 
mounting which could be adjusted vertically and laterally. His head was positioned 
so that the point mid-way between his eyes was in alignment with the axis formed by 
the Q apertures when swung into their observing positions. This could be judged 
from the foreshortened appearance seen by each eye when the two extreme Q 
apertures were in position simultaneously. The angular separation 6 between the 
Q axis and the four P apertures was 5°, 10°, 14° and 20° respectively. The P 
aperture selected for any particular experiment was illuminated by a light box and 
diffusing screen L,, placed immediately behind it; the diameter of the P apertures 
was 7/16 in. thus having the same angular subtense of 25 minutes as the Q aperture. 

The light sources in the two light boxes were 4-v. lamps whose intensity was 
controlled by a rheostat and voltmeter. The luminance of the apertures could be 
adjusted to suit the observer, but in most of the experiments it had a value of about 
1-5 candles/ft?. Flash exposure of the Q apertures could be arranged by intro- 
ducing a rotating sector at S. The sector was mounted on a gramophone motor 
with its axis on the same horizontal level as the Q apertures and the observer’s 
eyes ; this ensured that the light into the two eyes was intercepted simultaneously. 

It was considered that with the arrangements described, all monocular aids 
to distance perception other than accommodation had been eliminated. ‘To con- 
firm this, tests were made using one eye only (in this case the observing eye itself was 
brought into alignment with the Q axis), and no distance perception was found to 
exist over the range of 14 in. that could be studied with the apparatus. ‘This also 
confirmed that accommodation was a negligible factor in the observations. 


§3. METHOD OF OBSERVATION 


The task of the observer was to state whether, for the particular conditions of 
the experiment, the aperture Q introduced by the operator appeared in front (plus), 
or behind (minus), the comparison aperture P. ‘The reply was recorded by the 
operator in the column appropriate for that aperture, a new aperture was introduced 
and the observation repeated. A rapid sequence of observations could be main- 
tained as the apertures were swung in in random order until 25 replies had been 
recorded for each aperture used. Only a few of the apertures were required in any 
given experiment, the range and spacing to be covered depending on the distance 
sensitivity under the conditions of the experiment, as determined by a few trial 
observations. The percentage of plus and minus replies at each distance was then 
calculated and plotted against the distance. If the aperture Q appeared to be level 
with P, the observer was not allowed to state this, but had to make a guess of plus or 
minus. The level position was then, at least in theory, indicated by a 50%, 
plus, 50°% minus, estimate. 

Two conditions of observation, A and B, were employed. In condition A, the 
observer looked first at P and then at Q, so that both convergence and disparity 

T-2 
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contributed to his depth perception. The observer was asked to try to transfer 
his gaze from P to Q in one fixation jump, and certainly without any fixation pauses 
midway between P and Q. Eye movement studies (Lord and Wright 1949) 
suggest that in a quick movement, it is fairly easy to traverse most of the gap 
between two fixation targets in one jump, although a small residual correction may 
sometimes be necessary to complete the fixation on the new target. Ideally, the 
observer made his judgment of Q after a single switch from P to Q and he was 
definitely discouraged from looking backwards and forwards several times, 
otherwise the risk of a fixation pause between P and Q would have been too great. 

In condition B, P was fixated continuously and the position of Q had to be 
judged solely by extra-foveal vision. Convergence thus played no part and, since 
the disparity contributions to depth perception were the same for the A and B 
conditions provided the observer made no intermediate fixation pauses in the 
A observation, any contribution to the stereoscopic sense from convergence 
changes would be revealed by an enhancement of the stereoscopic acuity for 
condition A. The results for the various angles 6 enabled the relative contributions 
of convergence and disparity to be studied for increasing angular separations of 
the test objects. 

At 0=14°, however, a special situation arises since, when P is fixated, Q is 
imaged on the blind spot of the right eye, while when Q is fixated, P falls on the 
blind spot of the left eye. There is thus no possibility of retinal disparity nraking 
any contribution to the stereoscopic acuity, hence for condition A the acuity is 
due solely to the convergence factor, while for condition B, observation of Q 
reduces effectively to one-eyed vision. ‘This value of 0 was therefore chosen 
deliberately as a device to separate the convergence and disparity factors. 

Tests on the flash acuity were restricted to condition B. Aperture P was 
fixated and seen continuously, and the observer was required to judge the position 
of Q as seen by extra-foveal viewing during one, or at most two, flash exposures. 


§4. RESULTS AND DISCUSSION 


Typical results are shown in Figure 2, in which the percentage estimates are 
plotted against the position of the Q apertures for A and B observations at 6 =10°. 
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Figure 2. Stereoscopic acuity curves for 6=10°. Percentage estimates of location of aperture 
plotted against aperture position for conditions A and B. 
+ test aperture reported in front of reference aperture. 

— test aperture reported behind reference aperture. Observer: W. D. W. 


Curves as smooth as those in Figure 2, A, were quite common, but under 
conditions where the acuity was poor, as in Figure 2, B, the points were sometimes 


The Role of Convergence in Stereoscopic Vision 293 


more irregular. Some 150 or 200 observations were required for each curve, 
but the apparatus could be operated so quickly that these could be recorded in 
15 or 20 minutes. More extended observations, which might have been justified 
on statistical grounds, were undesirable visually if fatigue were to be avoided. — 
One method of specifying the stereoscopic acuity from the percentage curves 
would have been to give the displacement of the apertures necessary to pass from 
a 50% to, say, an 80% estimate. This might correspond to the displacement 
needed for the aperture Q to be just noticeably in front or behind P. A measure 
which is in some ways to be preferred is obtained from the slopes of the curves 
at the 50% level, as measured by the percentage change per inch, since this 
gives a very direct indication of the sensitivity of the observer to change in 
distance of Q. Using the latter measure, Figure 3 shows the variation of the 
acuity with 6 for both conditions A and B and for a number of repeat experiments. 
While the slope of the percentage curves should be regarded as measuring the 
acuity, it has been convenient in Figure 3 to use a logarithmic scale, since this 
gives a better idea of the spread of the observations. The spread is in fact fairly 
considerable, as must be expected from the nature of the observations. 
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Observer: W. D. W. 


At 6=14°, the slope was zero for condition B, since the observation was 
effectively monocular; this is indicated by a break in the curve. A discontinuity 
is also shown in the A curve because, while observations could still be made, 
the judgment had a different basis from that at the other three positions. 

Both the A and B curves show a decrease in the stereoscopic acuity with 
increase in 0. The decrease is more severe for condition B, no doubt because, 
for this condition, depth perception is more directly dependent on the resolving 
power of the peripheral retina. ; 

Some stereoscopic acuity could be recorded for condition A at @=14°, and it 
would seem that this must be attributed to the convergence factor. ‘The average 
slope value lies about 0-3 log units below the dotted portion of the curve, so that 
the stereoscopic acuity as measured by the slope is roughly halved by the removal 
of the disparity factor. It might then be deduced that, in the absence of the 
blind spot, convergence and disparity would make approximately equal 
contributions to the acuity at @=14°. 

This can also be seen from Figure 4, in which the ratio of slope A to slope B 
has been plotted against 6. Except at 9 =0°, slope A is always greater than slope B, 
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indicating that convergence is always making a significant contribution to the 
acuity and, since the ratio increases with 9, it would seem that convergence 
becomes steadily more important the greater the lateral separation of the objects, 
until at @=20° it has become the dominant factor. The extrapolation of the 
curve to show a slope ratio of 1:0 at @=0° is reasonable, since as the targets 
approach very close to each other, conditions A and B become virtually identical. 

The results shown so far have all been recorded by the author, who has 
repeated each set of observations some four or five times. Two other observers 
have also made a single complete series of experiments, giving results of the same 
general character, although differing in absolute sensitivity. Thus Figure 5 
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Figure 5. Stereoscopic acuity curves for 6=20°, for conditions A and B. Observer: R. N. W.- 


gives A and B percentage curves for 6=20° as recorded by R. N. W. and 
Figure 6 the A percentage curves for the blind spot condition, @=14°, as 
recorded by E. T. 

Flash acuity tests were made by W. D. W. for 6=5° for condition B only. 
The main purpose of the flash observations was to determine whether the small 
eye movements that may occur during apparently steady fixation make any 
important contribution to fusion and depth perception, on the lines suggested 
by Marshall and Talbot (1942). As shown in Figure 7, no significant 
deterioration in stereoscopic acuity occurred for the flash condition, and this 
result was confirmed for flash durations between 0-10sec. and 0-01 sec. “This 
conclusion is in apparent disagreement with Langlands’ results, but in the present 
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experiments the reference aperture P was being fixated all the time and this, 
no doubt, made the observer’s task easier. Nevertheless, the present arrangement 
appears to provide a legitimate test of the part played by small eye movements 
in the development of the retinal disparity contribution. The absence of any 
deterioration under flash conditions suggests that eye movements are not 
necessary to produce fusion of disparate retinal images and do not enhance 
stereoscopic perception. 

One effect experienced by W. D. W. during the blind spot experiment should 
be put on record. After a few minutes’ observation, the ability to locate the 
distance of Q relative to P sometimes disappeared almost completely, but after 
a minute or two’s rest, including looking round the room, depth perception 
was restored. Possibly some adaptation process occurs in the sensory fibres 
of the extra-ocular muscles, but the effect at least confirms that the depth 
perception recorded was not an artefact arising from fixation pauses intermediate 
between P and Q. 

§5. DEPTH PERCEPTIONMN DIPLOPIA 

Discussions on stereoscopic vision in terms of the retinal disparity theory 
generally seem to imply that fusion is an integral part of the depth perception. 
Yet experiments have been reported (Auerbach and v. Kries 1877, Tschermak and 
Hofer 1903) in which it was claimed that fusion is not required for the creation of 
the depth sensation ; a spark, for example, seen as a double image could be correctly 
located in front or behind a reference object. 

Under ordinary conditions of seeing it is almost impossible to test whether 
crossed diplopia can be distinguished from uncrossed diplopia without the sub- 
sidiary clues of parallax, perspective, etc., which nearly always exist, or without the 
assistance of momentary changes in convergence which affect the amount of crossed 
and uncrossed diplopia in opposite directions. Nevertheless, the possibility that 
the two types of diplopia can be distinguished in their own right, as it were,through 
their influence on depth perception, seemed so remarkable and so fundamental 
to an adequate theory of stereoscopic vision, that it seemed worth while to extend 
the present experiments to confirm this claim. 

|P 


jy. Observer 


1 


j4- 
“ ©) 


(a) Side Elevation 


Q, 


Qn, and Q,, Visual Axes : L 
(b) Plan 


Figure 8. Arrangement for testing depth perception in the presence of diplopia. 
(a) Side elevation. (6) Plan, showing location of double images. 


For this purpose, the apparatus illustrated in side elevation in Figure 8(a) was 
constructed. An illuminated reference aperture P was located at a distance of 
2 ft. 8 in. from the observer and 14 in. vertically above the line of the test apertures 
Q,andQ,. Q, and Q, were located at a distance of 2 ft. and 4 ft. respectively from 
the observer, but any one of five apertures of different sizes could be inserted at 
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either Q, or Q,. The angular subtense of P at the observer’s eye was about 
27 minutes of arc, while the Q apertures subtended angles ranging from 18 minutes 
to 36 minutes. No attempt was thus made to maintain a constant apparent size 
of the Q apertures ; instead, size was eliminated as a factor by having both small and 
large apertures at both the nearer and farther positions. 

In the experiment, P was fixated continuously, but Q, and Q, were only seen 
momentarily and for a single flash as the sector S rotated. This sector consisted of 
a cylindrical drum with a flat segment, which exposed the Q aperture when the 
segment crossed the line of sight. Suitable screening was provided so that, as 
before, the observer saw only the illuminated apertures in a black field. When the 
operator introduced either aperture Q, or Q, in random sequence both of size and 
distance, the observer had merely to indicate whether the aperture he saw in a 
single flash was the one in front or behind P. 

From the plan view of the arrangement in Figure 8(d), itis evident that when Q, 
is introduced, it will be seen as a double image projected at Q, (right eye) and Q;, 
(left eye) in the plane through P. When Q, is introduced, the corresponding double 
images will be projected at Qp, and Q;,. From the dimensions of the system, Qp, 
will be coincident with Q;,, and Q;, with Qg,, sothat the visual patterns are identical 
whether Q, or Q, is inserted. Since the apparent angle separating the double 
images was about 14° and they were viewed only about 2° below the fovea, fusion 
was manifestly impossible; the short duration flash (about 0-02 sec.) prevented 
convergence changes giving any clue as to whether Q, or Q, had been inserted ; the 
observer’s head was fixed, hence parallax effects were eliminated ; and no subsidiary 
aids from perspective, shadows, size, etc., were present. Yet the observer never 
had the slightest difficulty in deciding whether Q was in front or behind P and 
indeed the double image was seen clearly projected in front of, or lying behind, P 
with a true three-dimensional effect. 

Evidently, therefore, fusion is not an essential item in depth perception. 
While the observer cannot consciously tell which eye is seeing which of the double 
images (a matter of considerable interest and referred to by many writers, e.g. 
Helmholtz 1866, English translation 1925 p. 458) the information must neverthe- 
less be available to the brain and be utilized to give the position of the object in 
space. ‘The experiment supports the theory sometimes encountered in the litera- 
ture (see, for example, Carr 1935, p. 223) that depth perception, in part at least, 
results from a kind of mental surveying process in which Q, is located at the 
intersection of the lines RQ, and LQ,, (Figure 8(6)) mentally projected from the 
right and left eyes to their respective disparate images, and Q, at the intersection of 
RQpg, and LQ;,. The positions of Q, and Q, as judged in this way would clearly 
have to be related in absolute distance to the position of P as determined by the 
convergence angle between the visual axes PRand PL. The theory could then be 
integrated very satisfactorily with the idea of convergence as a significant factor in 
depth perception. 


§6. CONCLUSION 


Certain fairly definite conclusions emerge from these experiments : 

(1) The condition B observations indicate that retinal disparity, accompanied 
by fusion, contributes to a sense of depth. 

(1) The diplopia experiments indicate that retinal disparity, even in an extreme 
form unaccompanied by fusion, contributes to a sense of depth. 
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(ii) Whether fusion as such makes a contribution to stereopsis additional to 
that arising from the disparity itself cannot be decided from the present 
experiments. No doubt arguments for and against this possibility could be 
advanced. ; 

(iv) The flash experiments indicate that small eye movements are not a 
necessary item for fusion and do not enhance depth perception. 

(v) The blind spot experiments suggest that convergence by itself is an 
important factor in stereoscopic acuity. 

(vi) The relative values of stereoscopic acuity for conditions A and B at 
various separations of the targets suggest that convergence and disparity both 
contribute to the stereoscopic sense, and that the wider the field of view, the more 
important does the convergence factor become. 

Reference has been made at one or two points in this paper to the need to 
eliminate from the experiments subsidiary clues such as parallax, perspective, size, 
etc. This elimination was necessary to reduce the number of variables in the 
experiment, but there is ho suggestion that these aids are not genuine factors in 
three-dimensional vision. It may well be true that the sensation of depth was first 
acquired through binocular perception, but, either by association or otherwise, the 
monocular factors have come to be important items as well, and can by themselves 
give rise to three-dimensional seeing. Under ordinary conditions of viewing, 
stereoscopic vision is therefore the result of a mental analysis and synthesis of 
a quite complicated array of visual data. 
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Penetration by High-Velocity (‘Munroe’) Jets: I 


By D. C. PACK* ann W. M. EVANS 
Armament Research Establishment, Fort Halstead, Kent 


Communicated by N. F. Mott; MS. received 19th Fuly 1950, 
and in amended form 10th November 1950 


ABSTRACT. By means of certain simplifying assumptions a formula is developed for 
the penetration into a ductile target by a high-velocity (‘ Munroe ’) jet. The action of the 
jet is divided into two stages, each making its contribution to the total penetration. In the 
first stage a hole is formed by the lateral compression of the target as the jet penetrates it; 
the second stage begins when the last particle of the jet has ceased to act, the hole continuing 
to deepen until the residual energy in the target has been spent. At the high pressures 
set up by a Munroe jet the strength of the target plays only a subsidiary part in the 
phenomenon. 


§1. INTRODUCTION 


XPERIMENTS were carried out by Munroe (1888) in which gelatinous 
cartridges, each with a conical indentation in the end, were placed against 
a steel target and exploded. It was observed that these cartridges produced 
indentations in the target which were similar in form to those in the original 
cartridge. The discovery was not entirely new; at the end of the eighteenth 
century it had been known that the existence of a hollow in the forward end of a 
blasting charge led to an enhancement of the explosive effect along the axis. The 
‘hollow charge’ was used to some extent in mining in Europe during the early 
part of the nineteenth century, but it appears to have fallen into disuse, and 
Munroe’s later work led to little research of interest until recent years. 

Payman and Woodhead (1937) showed the importance of solid particles, 
carried in the detonation products, in producing and prolonging the intense 
‘end effects’. It is now known that these effects become even more remarkable 
when the indentation or hollow is lined with material ; a pronounced jet is formed, 
and the stream of material from the lining which acts as the penetrating agent has 
been called the ‘Munroe jet’. A considerable amount of pioneer research on 
phenomena associated with the detonation of an explosive charge with a lined 
hollow at the end away from the detonator was carried out by Evans and 
Ubbelohde (1950 a,b). A general summary of work on the formation of Munroe 
jets has been given in a paper by Birkhoff, MacDougall, Pugh and Taylor (1948). 
Experiments have provided valuable information concerning the behaviour of 
the Munroe jet, but many questions on the origin and manner of action of the 
jet remain unanswered. 

The purpose of the present paper is to investigate, from a theoretical point of 
view, the penetration by Munroe jets into solid targets. The mechanism of 
formation of the jet will not be discussed. It has been necessary to make some 
very sweeping simplifications and approximations in the formulation and solution 
of the problem, but it should be borne in mind that formulae which give only 
approximate answers may be of great assistance in practical work: it is well known 
by experimental workers that explosive charges made carefully to the same 
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specifications and from the same materials usually give results which exhibit 
an appreciable amount of scatter. 

Since the measurement of the actual physical conditions in a jet, even when 
possible, is a difficult experimental task, the object of the theoretical work has 
been to express results in such a form that the performance of a jet against 
different combinations of targets may be calculated from the values of certain 
parameters of the jet, these parameters being easily determined from the results 
of an initial small number of simple firings. This application will be discussed 
in a later paper. 


§2. DEPTH OF PENETRATION 


Evans and Ubbelohde (1950 b) measured the velocity of the jet from hollow 
charges lined with metal. The velocity of the head of the jet is of the order of 
10,000 feet per second, so that the pressure exerted by the jet, on any material 
which it penetrates, must normally be large compared with the plastic yield stress 
of the material. It follows that the velocity of penetration will be determined 
by the inertia of the target rather than by its plastic properties. On the other 
hand, the resultant deformation in the target will depend directly on the yield 
stress. The exact physical state of the metal in a Munroe jet is not yet known, but 
Evans and Ubbelohde speak of two types of jet, ‘fluid’ and ‘fragment’ jets, 
depending upon the nature of the metal (or other material) in the lining of the 
hollow. 

We assume, therefore, as the basis for a theory of penetration, that (a) the 
material in a Munroe jet is either liquid or solid for which the flow stress is 
negligible, or it is a jet of fragments which break up on impact; (0) the yield stress 
of the target material may be neglected, to a first approximation, in comparison 
with the pressures set up by the penetrating jet. 

It may be mentioned here (Evans and Ubbelohde 1950 a) that a metallurgical 
examination of steel blocks attacked by the Munroe jet from a hollow charge 
lined with steel showed that penetration was achieved by lateral expansion of the 
target material, the main displacement of the material being in a radial 
direction. Further, the diameter of the final crater produced in a target by a jet 
is larger than the diameter of the jet itself, indicating that the very high velocities 
of penetration lead to cavitation of the target material, as shown in Figure 1. 
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Suppose (Figure 1) that the material in a Munroe jet moves with uniform 
velocity V and that it makes a hole in the target of diameter greater than its own. 
The material of the jet, which breaks up on hitting the bottom of the hole, is 
pushed away; in the metallurgical examination referred to above, particles of 
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the steel lining were observed adhering in places to the walls of the cavity. 
Let U be the rate of increase of the length of the hole. Let the whole motion 
be in a steady state. If the jet has length L the time during which pressure is 
exerted on the bottom of the hole is L/(V—U). The depth p of the penetration 
is therefore given by 


paLu(Ve Cant: ly ang ikea iey aim (1) 


If O is the point in the axis of the jet at the bottom of the hole, U may be determined 
by equating the pressure of the jet at O to the pressure required to drive a hole 
with velocity U into the target. Axes are taken through O along and perpendicular 
to the axis of the jet; relative to these axes the jet moves with velocity V—U, 
and the material of the target with velocity U at infinity (see Figure 2). 

When U is very large, the material of the target may be treated as an inviscid 
fluid; the pressure that it exerts at O (the stagnation point) will then be $p,U?, 
assuming incompressibility ; p, is the density of the target. 

The pressure exerted at O by the jet will depend upon its structure. 

For a liquid jet or a solid jet of negligible flow stress the pressure will be 
$p,(V —U)?, where p, is the density of the jet. 

For a fragment. jet, at the high velocities observed, it cannot make any 
difference whether the fragments are of a hard or soft metal, liquid or solid; they 
will, in any case, shatter on whatever surface they hit and cause penetration by 
the momentum transferred. A fragment jet, for which the metal occupies only 
a part of the volume, will, under these conditions, exert a mean pressure equal 
to p,(V — U)*, p; being the mass per unit volume of the jet. 

Since for the opposite extreme (a fluid or solid jet where mass occupies 
the whole volume) the pressure at the stagnation point is half the above, it is 
reasonable to suppose that, for all jets of a given density p,, the mean pressure at 
the stagnation point will lie between the limits $o,(V — U)? and p;(V—U)?. The 
pressure may therefore be written as $Ap,(V — U)*, where 1 <A <2, or zp) (V— UP, 
where pi, =Ap;, the effective density: of the jet. Equating the pressure in the ye 
to that in the target at the stagnation point O, we have 


ppt =p) (VT ee (2) 
The velocity of penetration is thus 
U=Vjit+ fou = hee (3) 


For a fluid jet penetrating a target of the same material, the velocity of penetration 
is one-half of the velocity of the jet itself. 
The penetration is given by * 


p=LU(V—U)=L/(pilog. sas (4) 

With the chosen model, therefore, the interesting conclusion is reached that, 
to a first approximation, the depth of penetration by a jet into a given target is 
independent of the velocity of the jet, and depends only upon its length and density. 
For a given jet against different targets, the depth of penetration varies inversely 
as the square root of the density of the target attacked. 

It should be noted that these formulae apply only if the target is made of a 
substance of low compressibility ; for the very high pressures concerned, rubber 


* This result was published by R. Hill, N. F. Mott and D. C. Pack in a Ministry of Supply 
Report, January 1944. 
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or paraffin could not be treated as incompressible. Furthermore, the pressures 
exerted by the jet must be very high compared with the yield stress of the target 
material. This condition is usually fulfilled; if p, is put equal to the density of 


steel and V =10,000 ft/sec., 4p,V? = 2,400 tons/in®, while the yield stress of steel 
is of the order 40 tons/in®. 


§3. LIMITATIONS OF PRECEDING THEORY 


There are three important cases for which the simple theory outlined above is 
insufficient to account fully for the penetration observed. 

From (4) it would be expected that penetration into armour would be the same 
as into steel, since their densities are equal, but it is an experimental fact that the 
depth achieved is definitely less. The difference must arise from strength 
considerations; in obtaining (4) the yield stress of the target was neglected in 
comparison with the pressures set up by the high-velocity jet. 

The second difficulty arises when penetration into lead is measured. By 
equation (4) the penetration by a given jet is inversely proportional to the square 
root of the density of the target, and should therefore be less for lead than for 
steel; experiments on lead, however, have resulted in penetrations approximately 
double those which would have been predicted on theoretical grounds. 

Lastly, the formula gives a value for the penetration which is independent of 
the distance between the base of the cavity charge and the face of the target 
(the “stand-off’), while it is known from experiment that there is a marked 
dependence upon this distance. The effect of stand-off is due principally to a 
velocity gradient which exists in a Munroe jet, by virtue of which the jet elongates 
as it travels. The authors have shown however (in an unpublished paper) that 
if the assumption be made that there is a linear velocity gradient in the jet, then 
for targets of steel, armour, lead and aluminium, provided the stand-off is not 
too small, results taken at a given stand-off may be analysed sufficiently accurately 
by means of the simple formulae developed in the present paper. 


§4. THE EFFECT OF TARGET STRENGTH 


To obtain the fundamental equation (4) it is assumed that the yield stress of 
the target may be neglected in comparison with the pressure developed by a 
Munroe jet during penetration; that is, Y/p,V* is neglected in comparison with 
unity, where Y denotes the dynamic yield stress of the target and V is the velocity 
of the jet. To make a correction for the effect of strength, one may suppose 
that the penetration is expanded in a series of powers of the non-dimensional 
parameter Y/p,V?, the first approximation being 


P=L(p; lpr) Pa Vip) tee (5) 


where «, is a function of the densities of the jet and the target. 

Estimates made by Evans and Pack (1951) show that for armour, which has 
_a very high yield stress, the correction term a, Y/p;V* can be as much as 0-3. ‘The 
effect of strength in such a case is to reduce the penetration given by (4) by 
30 per cent. In lead the correction term is negligible for the same jets. ‘The 
adjustment in the ratios of penetration in lead and armour caused by the strength 
correction is however, not sufficiently large to account for the relative excess of 
penetration observed in lead. 
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§5. SECONDARY PENETRATION 


Secondary penetration is defined as that which occurs after the main jet 
penetration is ended. 

It is known experimentally that a Munroe jet makes a wide hole in a target 
of lead, and this is attributable to the low yield stress of this metal. The jet sets 
the lead in motion laterally, and the final width of the hole represents the diameter 
at which the energy of the motion has been expended against the resistance of the 
target to the flow. This being so, we may consider what is the position in the 
target when the last particle of the jet has been consumed. A high pressure has 
been applied to the bottom of the hole, and suddenly ceases to act. The hole 
continues to expand until the motion begun by the pressure has been damped 
out by the resistance to the flow. If the previous pressure be uniformly 
distributed over the bottom of the hole, the latter, expanding freely, adopts a 
hemispherical shape. This has been found to be roughly the shape of the bottom 
of the crater made by Munroe jets in general, when plug effects are not present 
(Evans and Pack 1951), so that the increment in penetration due to the free 
afterflow of the target is approximately equal to the radius of the hole made by 
the jet. 

Under these conditions the penetration of a Munroe jet into a ductile target 
may now be written in a form which allows for the effect of strength and for the 
secondary penetration. Denoting by p the total penetration into a target and 
defining the elementary penetration p’ as the quantity given by the equation 


p =Lyv/(p; /px)s then p=p'(1—«,Y/p,V?) +7, Sec (6) 


where r is the secondary penetration, which will be put equal to the radius of the 
hole made by the jet. The quantity p’(1—«,Y/p,V) will be called the primary 
penetration. 

In another paper it will be shown how the two phases of the penetration 
may be separated experimentally, and how the formuJae developed above may be 
used for practical purposes. 
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ABSTRACT. Work is described in which the two stages of penetration by a high-velocity 
(‘Munroe’) jet were separated experimentally. The large penetrations measured in lead 
targets are shown to result from the flow which takes place in the metal after the jet itseif 
has been consumed. For a given jet at a given stand-off it is possible to predict the 
penetration into a combination of targets from the results of a very small number of standard 
experiments. The method depends upon the calculation of a certain quantity which is 
constant for a given jet at a given stand-off, and examples are given of the determination of 
this quantity from experimental data. 


§1. INTRODUCTION 


N a previous paper (Pack and Evans 1951, to be referred to as I) is given a 
] theory of penetration by Munroe jets in which the penetration observed in 
a ductile target is divided into two parts: (i) the primary penetration, which 
depends chiefly upon the density of the target, and (ii) the secondary penetration, 
which takes place after the jet has ceased to act, and which represents the 
expenditure of the energy left in the target at the time the last particle of the jet 
has been consumed. In the present paper work is described in which the factors 
involved have been separated experimentally, and an analysis of the results provides 
reasonable confirmation of the proposed theory, making it possible to predict the 
penetration of a high-velocity jet into any combination of ductile targets, using 
the results obtained from a very small number of standard experiments. 

It is known that when a conical metal lining is used for a hollow charge, the 
metal divides into two parts, the high-velocity Munroe jet in front and a slower, 
much wider (and often solid) piece of metal, usually called the ‘plug’, behind. 
The plug is capable of increasing the penetration, acting after the jet has been 
consumed, and our experiments have thrown some light upon its behaviour. 


§2. THE RELATION BETWEEN ENERGY AND VOLUME DAMAGE 

It has already been mentioned (I) that a Munroe jet displaces the target 
material almost laterally during penetration. The work done is therefore 
approximately equal to that required to make a cylindrical hole. It has been 
shown by Bishop, Hill and Mott (1945) that, for static punching, R, the work 
per unit volume required to form a hole deep in an infinite block, is, to a good 
approximation, independent of the shape of the hole. For materials which do 
net harden appreciably R is of the order of 3:5 to 4 times Y, where Y is the initial 
yield stress. For annealed metals R may be between 5Y and 6Y. When a hole is 
made dynamically R depends also upon the mean rate of strain. For large strain 
rates, however, R will not vary between wide limits of the velocity of penetration, 
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on account of the logarithmic dependence of yield stress upon rate of strain. 
‘Thus, neglecting surface effects, 


work done per unit length of hole = R x area of cross section of hole ...(1) 


This relationship is presumably also valid for the dynamic penetration by a 
Munroe jet, as long as we are sufficiently far from the surface of the target. The 
latter restriction is not very important for Munroe jet penetrations, since they 
are usually deep compared with the radius of the surface crater. The equation 
therefore provides a connection between the energy given up by the jet and 
the volume damage in a target. 

The equation for the penetration due to a cylindrical jet of mean density p; 
and length L entering a target of density p, and yield stress Y was given as 

p= LApy|px)#2(1 — 04, ¥ /p,V2) + = L(Apylon)!2(1 —RR) +7, using (1), ...(2’) 
where A is a factor between 1 and 2 depending upon the nature of the jet, k is a 
factor depending upon the properties of the jet and the density of the target, 
and 7 is the radius of the crater (measured near the bottom of the hole). r is the 
secondary penetration; L(p,'/p,)"?, where p;’ =Ap;, is the elementary penetration 
obtained by neglecting both target strength and secondary effects. 

For the sake of abbreviation, (2’) is written 
p=p MRR Vi ie eee (2) 
p'(1—RR) is called the primary penetration. 

The ratio of the values of R for two targets may be found, by virtue of (1), by 
comparing the ratio of the volumes of the craters made by a standard jet. For 
steel and armour, which, having equal density, have the same factor k in (2) and 
the same elementary penetration, the effect of the strength on the penetration may 
therefore be calculated, without the assumption of any numerical details 
concerning the jet, by solving two equations of the type (2) for the unknowns 
p’ and kR, the equations being supplied by experimental observation of p, 7 and 
the volumes of the holes made by a given jet. . 


§3. PREVENTION OF AFTER-FLOW: EXPERIMENTAL WORK 

It follows from what has already been said that the penetration would be 
wholly primary if ‘after-flow’ or secondary penetration could be prevented. 
The most simple test of this theory is to back the target in which secondary 
penetration is considered to be of importance with a second metal, such as steel 
or armour, for which the resistance to flow is very high. We may then analyse 
the results of the combined penetrations and compare them with those observed 
in massive blocks of the individual metals. 

A series of experiments on combined targets has been carried out to examine 
these hypotheses. Charges were made up consisting of 60 gm. of explosive 
confined in 1g in. internal diameter mild steel tubing of + in. wall thickness, 
and lined at the base by (i) 80° pressed mild steel cones 0-032 D thick, and 
(i1) pressed caps of tin of 0-73 D radius of curvature and 0-044 D thick, where D is 
the internal diameter of the charge. The targets used consisted of massive blocks 
of mild steel," homogeneous armour plate and lead, various thicknesses of lead 
backed with mild steel, and various combinations of armour and mild steel. For 
the conical type of charge the stand-off (i.e. the distance at which the base of the 
charge is placed from the target) was D, and for the tin caps 2:3 D. 

The results are given in Tables 1 to 5. 
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Table 1 
Mean Mean volume Mean radius 
penetration... of hole of bottom 
(mm.) (cm?) of crater (mm.) 
Mild steel cones 
(allowing for effect of plug) 
Massive mild steel 66:5 11 8:5 
Massive homogeneous \ 61:5 7 6 
armour 
Massive lead 83 382 37°5 
Mild steel laminae 
stacked together \ oe ue) $3 
Tin caps 
Massive mild steel 65 10 5 
Massive homogeneous \ 55 6 3 
armour 
Massive lead 97 472°5 40 


Table 2. Mild Steel Cones 


Thickness of Penetration 
lead penetrated of steel underneath K TE 
(mm.) (mm.) 
Lead backed with mild steel 10 49 155 180 
plates 10 49 1155 180 
15 35 130 155 
O77, 18 153 179 
37 25 . 174 200 
39 19 163 188 
49 7 165t 186 
54:5 2 Ey 
60 2 * 
62:5 0 : 
Stack of mild steel plates 172 197 
Massive lead 154 280 
Mean 158 
Standard deviation 12 


+ Assuming zero primary penetration into steel. 
* Values omitted for reasons given on page 309. 


Table 3. Mild Steel Cones 


Thickness of Penetration of 
| armour penetrated steel underneath K IR 
(mm.) (mm.) 

Armour backed with mild 19°5 37 145 iil 

steel plates 19-5 40 154 180 

39°5 DP) 163 188 

39-5 25 172 197 

39°5 20 157 182 

Stack of mild steel plates 7 197 

Massive armour iP 190 
Mean 162 
Standard deviation 10 

U 
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Table 4. ‘Tin Caps 


Thickness of Penetration of 
lead penetrated — steel underneath K TG 
(mm.) (mm.) 
Lead backed with mild steel 14 44 177 194 
block 27 38 201 218 
38 27 201 218 
51 12 + 195 212 
56 6:5 194 211 
75 @) “ 

Massive steel 199 216 
Massive lead 2192 328 

Mean 194 

Standard deviation 8 


* Value omitted for reasons given on page 310. 


Table 5. ‘Tin Caps 


Thickness of mild Penetration of 
steel penetrated armour underneath IK 1s 
(mm.) (mm.) 
Mild steel backed with 35 26:5 206 218 
armour block 39 25 214 225 
50 12:5 203 214 
62 0 = 
63 0 = 
Massive armour 199 210 
Massive steel 199 216 
Mean 204 
Standard deviation 6 


* Values omitted for reasons given on page 310. 


§4. PENETRATION PHENOMENA : 

It will be of interest at this point to remark on some of the observations made 
during this series of experiments, which bring to light characteristics of the 
behaviour of the plug from a conical lining, the plug being the more slowly 
moving but larger mass of metal which follows the high-velocity jet. No plug is 
found when caps of tin are used. 

It has been observed, by means of x-ray flash photographs, that the plug 
from a conical lining is much wider than the jet. In consequence, it often happens 
that, with a steel target, in which the jet bores a hole of only 4—5 times its own radius, 
the plug becomes trapped in the hole, particularly if it is not well aligned with 
the jet. In penetration of lead, however, the width of the hole is so much greater 
that the plug is enabled to travel freely until the jet has ceased to act. The plug 
is then to be considered as a projectile, and its penetration has been measured in 
the way outlined below. 

For a massive lead target, attacked by the jet from a conical steel lining, the 
shape of the bottom of the crater was found by partially filling the hole with 
paraffin wax in which a T-shaped handle was inserted. On cooling of the wax, 
the cast was removed by pulling the handle. A photograph of the cast is given in 
Figure 1(b)*. By rounding off the contour at the discontinuity, the extent of the 


* See Plate at end of issue. 
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total jet penetration could be judged, and it was then found that the plug ~ 
penetration was approximately 30mm. For a lead block of critical thickness 
(i.e. of thickness equal to the total penetration for lead alone), backed with steel, 
it was observed that there was a sudden constriction of the hole at the bottom of 
the lead block, together with a penetration of 5mm. into the steel. The latter 
measures the penetration of the plug into steel, and due allowance for this factor 
has been made in the calculations on cones. 

The resistance offered by a target to the passage of a projectile is dependent 
to a high degree on the shape of the projectile and upon its deformability; for 
the plug we may perhaps prescribe limits by consideration of the laws found by 
Hill (1944) for resistance to ogival-shaped heads. For a given projectile there 
is a critical velocity above which cavitation sets in, ie. the target ceases to 
maintain contact with the projectile along the whole of its length. Below the 
critical velocity resistance is proportional to Y, the dynamic yield stress of the 
target, and an easy calculation shows that penetration is inversely proportional 
to Y. At velocities which are very high compared with the critical velocity, the 
resistance is proportional to pV”, where p is the target density and V is the velocity. 
Under these conditions penetration is inversely proportional to density. 
Comparing steel and lead, we see that we should expect the penetrations by the 
plug to lie between near-equality (corresponding to the density law) and an 
excess in lead of proportion equal to the ratio of dynamic yield stresses of steel 
and lead; in fact, the observed penetration of.the plug in lead is about six times 
what is observed in steel, so that the ratio lies between the limits prescribed. 

The contour of the hole in a massive lead block caused by a tin cap is given 
(Figure 1(a)) for the sake of comparison with that caused by a coned charge. 


§5. SURFACE CHARACTERISTICS IN A COMBINED 
TARGET ATTACKED BY A PROJECTILE 

Interesting as further experimental evidence on the behaviour of a plug as a 
projectile is the fact that for coned charges against critical thickness of lead 
resting on steel (or more generally where the thickness of lead lies between this 
value and—a lesser value—that of the primary penetration) the steel shows 
an upward lipping. ‘The latter is well known as a characteristic of projectile 
penetration. When the jet itself succeeds in penetrating into the steel backing, 
the front face characteristic in the steel is a depression. 

Dr. G. O. Baines, of the Armament Research Department, kindly verified 
that the upward lipping of the lower plate occurs in projectile penetration when 
lead rests upon steel, just as for plug penetration. He fired bullets of hard steel 
with 4 calibres radius head, coated with copper on the body, with velocities 1,500 
and 1,900 ft/sec. into lead blocks, 1}in. thick, standing upon steel plates. Two 
different hardnesses of steel plate were used. As far as could be seen, there was 
no difference in the front surface of the steel plates with or without lead, due 
allowance being made for the loss of velocity in the lead and the skew penetration . 
of the bullet after passing through it. The lead was found to have been pressed 
hard against the steel plates AA (Figure 2), and a cavity BB was formed which 
accommodated the front lip of the steel. 

It is seen that, as the head of the projectile penetrates the steel, the tendency to 
lip succeeds on account of the pressure of the steel against the softer lead, 


maintained by the presence of the projectile. We know, also, that lipping occurs 
u-2 
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if the lead is moved away by cavitation before the steel is attached (plug 
penetration). It appears, then, that for jet penetration the velocity of free 
expansion of steel must be greater than that of lead, or at least of the same order, 
and that the motion of the steel is then inclined to be in the direction of motion 
of the jet, as is observed in the interior of a block. In other words, the jet does 
not meet a free steel surface. 


§6. MULTIPLE TARGETS 


The laws which hold when a Munroe jet penetrates a series of targets of 
varying densities and strengths may be derived as follows: 
Suppose the jet to penetrate completely (n—1) targets of thicknesses py, 
Po +++» Pnuy densities p,, po, - - +» Pn—1, the strengths of which involve terms 
R,Ry, RoRy . . . » RyRy in equation (2). Let the jet penetrate only partially, 
a distance p,, into the mth target, of density p, and strength correction k,R,,. 
Let the secondary penetration in the last target be 7,, and suppose that the targets 

Direction 

of Bullet 

Front. of Lead Block 


A B B A 
Steel Backing 


Figure 2. 


are so arranged that there will be no secondary penetration in the first (n—1) 
targets. Then if successive lengths L,, L.,..., L, of the jet are used up in 
penetrating successive targets, the whole jet length being L, we have the following 
series of equations: 


Pi=Ly(p; /p1)"7(1 — RR), 
D2 = Lo(p;'/p2)"(1 — keRs), 

Praz Ln-1(p; |Pn—1)"71 = Ng p tng) 
Pn= L,(p; |p»)! = RRs) +1 yn) 


from which it is easy to see that 


.: Pin V Pm ’ Y “/ p 
>3 Ci ¥ n nN 
m=1 1 —ki Rin Ly/p; 1 1 = k,R» ye ee MOtichon dd (3) 
S  Pin/ Pm TaN/ Ba 
Se) ech A Om SS v Sem 
or ail R Rin aceree Kyo & (0 Sige (4) 


where K is an absolute constant for a given jet at given stand-off. 
For the double targets used in the experiments of this paper we use the form 


Piv Pi ri PoV Pa Va Pa aK 
| 1—k,R, 1~k,R, 1-&,R, 
and evaluate K in order to assess the reliability of the theory. The values obtained 
(after allowance for the plug in the case of the steel cones) are included in the 
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Tables, and cover the range from a semi-infinite block of the front material to 
zero thickness. The effect of the removal of the backing is clearly seen in the 
values of the expression 


Rie PivP1 PoV Po 
TSR oO TORR’ 
which would be constant in the absence of secondary penetration. The leap 
in K’ accompanying the transition from a block of a certain limiting thickness in 


lead backed by steel to a massive block of lead is striking evidence of the 
secondary effect (Tables 2, 4). 


§7. EXPERIMENTS WITH MILD STEEL CONES 

The method of application of the formulae may best be seen by means of an 
example. 

The appropriate values of kR are found from the results on massive targets. 
It is seen, for example (Table 1), that the jet from the mild steel cone penetrates 
66 mm. of mild steel plates, with crater volume 10-5 cm? and mean radius 8-5 mm. ; 
in homogeneous armour plate the respective values are 61:-5mm., 7-0.cmi, 
60mm. Since the densities are the same, k is the same for each. The ratio 
of Ry (steel) to Ry, (armour) is inversely proportional to the volumes; thus 
R,/Rz=1-5. Hence we have, using (2), 


66 =p,'(1—RgRg) + 8-5; 61:5 =py'(1—-1-5kgRy) +6-0 
with p,’=p,’, on account of the equality of the densities of armour and steel, 
from which we deduce that p,’ =p,’ =61-5mm., kgR,=0-065. Thus, in the mild 
steel plates for example, we see that the elementary penetration is 61:5 mm., the 
primary penetration is 57-5 mm. and the secondary penetration is 8-5 mm., the 
last two combined giving the total observed penetration of 66 mm. 
The strength of lead is so small compared with that of mild steel that the 


term kR may be neglected in this case. Hence the expected penetration into a 
massive lead block will be 


P=Ps (Pstect/Preaa)”? + radius of crater in lead =50-9 + 37-5 =88-4 mm, 


The observed mean penetration was 83 mm. (Table 1). 

It is to be noted from the Tables that the mean value of K is practically the 
same for the two sets of experiments, and that the standard deviation is 8% and 
6% of the mean respectively. Certain values are omitted in the computation of 
the means, since they correspond to thicknesses of lead greater than the primary 
penetration; that is, there is some secondary penetration in the lead in these 
cases. The standard deviation lies well within the experimental consistency 
obtainable; it is possible to find as much as 20% variation in the measure of a 
penetration between two identical rounds under identical conditions. The 
tables of K can, then, be taken as satisfactory evidence in favour of the theory. 


§8. EXPERIMENTS WITH TIN CAPS 
The reason for the use of the tin caps was to obtain a set of results which 
would be free from the liability of error introduced by the presence of a plug. 
It has already been remarked that the contours of the holes made by these jets 
had no discontinuities ; the results of the experiments were therefore much easier 
to interpret correctly. Furthermore, it was possible to use blocks of the backing 
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material instead of plates, as in the previous case, since it was no longer necessary 
to check the true position of the plug; this avoided errors due to discontinuities 
which might occur in a stack of plates, though from Table 1 this error is seen to 
be small. 

Denoting the effective resistances to flow of steel and armour by Rg, Ry, as 
before, we have R,/Rg~1-7, ps’ =71-4mm., kgRg=0-16, RyRy =0-27. We see, 
then, that the primary penetration of the tin cap was 84%, of the elementary 
penetration into mild steel, and 73° of it into armour. 

The corresponding penetration expected theoretically in lead is 99 mm., and 
the mean observed penetration (Table 1) was 97mm. The mean value of K is 
somewhat higher for the steel-armour combination than for the lead-steel, but 
again the difference and the standard deviations are small. 

As before, certain results have been omitted in the calculations since the 


target thickness of the front material exceeded the primary penetration for the © 


material. It is of interest to observe here, as confirmation of the hypothesis 
that afterflow occurs after the jet process has ceased, and is not a simultaneous 
phenomenon, that for the above targets, although the upper material was 
completely penetrated, there was no effect on the backing material—this, too, 
despite the fact that the thickness of the front target was not as great as the total 
observed penetration in a massive block of the same material. 
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ABSTRACT. In order to extend the study and use of the elastic properties of geological 
bodies to include effects involving the amplitude and frequency of elastic radiation, a new 
source of energy is needed. 

A large moving-coil transducer is described which, when attached rigidly to the rock, 
radiates sinusoidal vibrations of controllable amplitude, frequency and duration. The 
operation of the vibrator is analysed in detail by the method of electromechanical analogy. 

For experiments in a chalk-mine, energy at each of the frequencies 300, 600 and 1,000 c/s. 
was radiated in the form of a square-enveloped pulse, of 20 msec. duration. Vibrations were 
detected in the rock by means of a piezoelectric microphone, amplifiers and an oscillograph, 
with an overall magnification of about 10°. Oscillograms show the complex disturbances 
produced by the arrival of direct and reflected energy, at distances out to the limit of the mine 
at 450 ft. Measurements of amplitude at 130 ft. are used to examine the mechanism of the 
transfer of energy from the vibrator to the rock. 

The power radiated was about 0-05 watt at 600 c/s., of which about one-sixth was in the 
compressional mode and five-sixths in the shear mode; but this power may readily be 
increased. 


Si INT RODUC TLON 
UR present knowledge of the phenomena involved in the propagation of 
elastic energy in the ground has been gained from the study of natural 
and artificial earthquakes, and is limited to such aspects as are révealed 
by energy in the form produced by these catastrophic sources. Seismic and 
explosive waves are characterized by a high intensity, rapidly established, 
followed by a complex of vibrations of indeterminate form. From measurements 
of the time at which the impulse, defined by its onset, reaches known points of 
observation its path may be calculated, and hence the geological structures 
traversed may be accurately mapped. ‘Thus seismology explores the earth’s 
deep interior, and seismic prospecting the outermost few miles of the crust. 
Since no such precise information is conveyed by other features of the 
impulse, phenomena involving the frequency and intensity of the propagated 
energy have been used only qualitatively, or studied as disturbing effects. 
Sufficient concerning these phenomena has already been discovered, however, 
to encourage a special study of them. For such a study the frequency content 
of the radiated energy must be brought under control and much simplified. 
Efforts to control the frequency of explosive waves have had little success; an 
essentially more refined source is clearly required. At the same time, to dispense 
with the well-defined onset of the seismic impulse would be to abandon the 
means of determining the path followed, upon the knowledge of which the 
study of many aspects of propagation depends. For some purposes, on the 
other hand, continuous vibrations are useful. The apparatus described here 
radiates sinusoidal vibrations at various audio frequencies, either continuously 
or in the form of a square-enveloped pulse. 
* Now at Geophysics Department, D.S.I.R., Wellington, New Zealand. 
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§2. PREVIOUS EXPERIMENTS AND APPARATUS 


“The most successful published work on the radiation of elastic energy into 
the ground from a mechanical source was that of Howell, Kean and Thompson 
(1940). A moving-coil apparatus, weighing some 300 lb., was used to generate 
vibrations of frequency 400 c/s., which were detected out to a distance of 1,200 ft. 
at the surface and down to 100 ft. ina borehole. More general interest is attached 
to the measurements in depth than to the surface waves. ‘The borehole 
measurements were extended to a depth of 400ft. with a magnetostriction tube 
source, set into resonance at 900c/s. by means of a condenser discharge. The 
experiments were carried out below the water-table in ground consisting mainly 
of water sands (Howell, private communication), and since water also provided 
the connection between the vibrators and the ground, the medium probably 
resembled water itself rather than solid rock. 

In experiments on a smaller scale two different methods of coupling a vibrator 
rigidly to rock have been used. Haskell (private communication) achieved a 
range of 165 ft. in limestone, using a small magnetostriction rod furnished with a 
tapered brass end, which fitted tightly into a cylindrical hole in the rock. The 
magnetostriction rods used by Bradfield (1947), on the other hand, had flat ends, 
to give contact with the medium over a circular surface, but it was found more 
convenient to cement suitably designed metal seatings into the rock than to grind 
the surface of the rock itself. Both these workers were concerned with frequencies 
of several kilocycles per second. 

At very low frequencies vibrating devices, such as the eccentric flywheel, 
have been widely used to test the ground as a foundation for engineering structures. 
A measurement is usually made of the frequency of maximum ground 
displacement ; this is often called the ‘ natural frequency of the ground ’ (Crockett 
and Hammond 1949), although it is determined by the stiffness of the ground 
together with the mass of the machine. In general, the propagation of energy 
is of no significance to the measurements. In some cases, however, a kind of 
resonance may arise in a layer of rock of appropriate thickness, or perhaps from 
a complex set of reflected waves arriving from various depths and directions 
(Byerly, private communication). 

Sharpe (private communication) has attempted to study the transmission 


properties of surface rocks, at frequencies of 10 to 200c/s., by means of a 


variable-reluctance steady-state vibrator. Severe attenuation and the difficulty 
of coupling the detector to the rock limited the range of measurement to 12 ft. 
In seeking a method for measuring the thickness of overburden, Northwood 
(private communication) tried various electromechanical vibrators, the largest 
of which was an array of rochelle-salt units, about 1 ft. in diameter. For this 
problem, however, blasting caps proved satisfactory, though success was confined 
to waterlogged ground. 

Little progress has evidently been made towards a general study of elastic 
propagation by means of electromechanical vibrators. In the work described 
above, special difficulties have usually arisen out of the particular objects in view. 
Surface rocks, though of obvious importance to practical methods of prospecting, 
are notoriously poor elastic media, and work at the surface has been diverted from 
the main task by the convenient use of water. A preference for portable apparatus 
has tended to limit the power of sources or to encourage the use of excessively 
high frequencies. ‘Thus the major preliminary problems have remained: to 
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generate a large mechanical energy at suitable frequencies and to impart a useful 
fraction of this energy to rock so that it will be radiated as elastic waves. 


§3. THE VIBRATOR 

In order to radiate elastic energy in rock, by contrast with fluid media, a 
relatively large force must be developed, since the displacement is necessarily 
very small. ‘The frequency of vibration for the present experiments was chosen 
to give a wavelength of the order of ten feet, so that discontinuities at ranges of a 
few hundred feet could be adequately resolved without incurring too great 
attenuation. ‘Thus a reasonably transportable apparatus was required, capable 
of producing large pressures at frequencies around 500 c/s. A convenient nucleus 
for such an apparatus was provided by an Admiralty Fessenden oscillator, of the 
type formerly in service for underwater signalling and depth-sounding. This 
device was suggested by Dr. A. B. Wood, and a machine was kindly made available 
by Messrs. de Havilland Propellers, Ltd. 

The Fessenden oscillator (Figure 1) is a transducer of the moving-coil type, 
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Figure 1. The vibrator. 


with the special feature, of great advantage for the present woik, that the moving 
part is a large, robust cylinder of copper. The cylinder is suspended in the 
air-gap of an annular electromagnet, on the core of which is wound an a.c. coil 
only slightly smaller than the copper cylinder. Close inductive coupling is thus 
achieved between the a.c. coil and the cylinder, in which large currents may be 
induced, and, since the cylinder lies in the pD.c. field, a corresponding mechanical 
force is set up between the cylinder and the electromagnet. 

For use in water the oscillator had incorporated a massive diaphragm, some 
30in. in diameter. With this diaphragm, and its mounting-flange removed, the 
bulk of the oscillator was reduced virtually to that of the electromagnet alone, 
consisting of a steel cylinder 8in. high, 20in. in diameter and weighing some 
6001b. The weight of the enclosed copper cylinder was about 61b. The. 
mechanical problem was to connect these masses to each other and to the rock 
so that a considerable part of their energy of motion would be communicated to 
the ground and radiated. 

The analysis given below leads to three practical conclusions: first, the copper 
cylinder should be rigid with the rock; secondly, the magnet mass should be 
supported by springs of suitable compliance; thirdly, the force should be 
distributed over as wide an area of the rock as practicable. The first and third 
requirements are somewhat contradictory, and called for a compromise in the 
design. ‘To specify the best conditions with any accuracy would involve the 


314 F. F. Evison 


calculation of every effective compliance present, an excessively complicated 
task in a system of such high mechanical impedance. The simple precaution was 
therefore taken of using massive steel members wherever rigidity was desired, 
while the diameter of the coupling system was made only double that of the 
copper cylinder. Measurements will be quoted later to show that the coupling 
system, constructed in this way, possessed adequate rigidity and effected a 
satisfactory transfer of energy. 

The magnet mass was mounted, by means of three springs, upon a circular 
steel plate, 1}in. thick and 18in. in diameter. In order to leave room between 
magnet and mounting plate for the springs, and also for the power leads, a steel 
ring of length 3in. and wall-thickness 3in. was inserted between the lower end 
of the copper cylinder, which protruded slightly from the ends of the magnet, 
and the mounting plate. This ring fitted tightly to the mounting plate, and 
accommodated the copper cylinder in a shallow groove in its upper end. The 
copper cylinder could then be clamped between the mounting plate below and 
a $-in.-thick steel disc above, by means of a steel rod of diameter 1}in., fitted 
with a nut at each end, passing through a hole down the centre of the iron core of 
the magnet. In this way the copper cylinder was fixed as securely as possible to 
the mounting plate, to which both upward and downward forces could be 
transmitted. 

To anchor the mounting plate to the rock, three massive steel bolts were 
fitted at equal intervals around the edge. The bolts, which were 2} in. in diameter, 
were designed for cementing into holes in the rock 3in. deep. With the holes 
tapered upwards, and using a fluid mix of two parts of washed sand to one of 
cement, a good connection was achieved, which remained sound throughout the 
experiments. ‘The essential mechanical features of the vibrator may be seen 
in Figure 1. 

Electrically the Fessenden oscillator consists of an iron-cored transformer 
with a single turn for secondary, together with a p.c. magnet. The A.c. was 
derived from an alternator driven by an electric motor. Both the motor and the 
D.c. supply were operated from 230V. a.c. mains. 


§4. MECHANICAL ANALYSIS 

Theoretical formulae have been calculated for the radiation impedance, in a 
solid medium, of various ideal types of small radiator (Bradfield 1947). Since the 
present experiments were carried out in the narrow galleries of underground 
quarries, where surface waves could not arise, the vibrator may be regarded as 
an oscillating sphere immersed in the rock, but coupled over the lower half of its 
surface only. Halving each term of the formula, quoted by Bradfield, for a 
completely immersed sphere, the radiation impedance of the vibrator becomes 


Z, = 8n9r4f2(2/V,3 —1/V,2)VV,d/3 +j(20r3fd/3—4rV.V.dlf), .....- (1) 


where r is the effective radius of the vibrator, f is the frequency, V, is the velocity 
of shear waves in the rock, V,, that of compressional waves, and d is the density 
of the rock. Here the positive reactive term clearly involves the mass of rock 
set bodily into motion, the negative reactive term indicates the stiffness of the 
rock, and the resistive term represents the ability of the rock to dissipate energy 
by radiation. Expressed in terms of mechanical elements, the formula may be 


itt 
eta Z,=1/G,+jloM,—1eaC.), ee) 
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where G,, M,, C, are respectively the effective leakance, mass and compliance 
of the rock and w=2zf. Some opportunity is given by measurements on the 
propagated energy, as described later, to test the validity of formula (1) in the 
present work. 

A mechanical analysis of the vibrator has been carried out, with the help of 
electrical circuit theory, by means of a formal analogy between mechanical force 
and electric current and between relative velocity and electric potential 
difference. This analogy, which compares electrical impedance with mechanical 
admittance, defined as the ratio of velocity to force, is preferable to the more 
usual analogy comparing impedances in both cases for two reasons. First, the 
mechanical system is much more readily represented as an admittance diagram 
than an impedance diagram (Gehlshoj 1947); secondly, the analogy conforms 
with the fundamental electromechanical relations for a moving conductor. 

The stages in the analysis are shown in order in Figure 2. Figure 2 (a) 


(Cc) d) 


Figure 2. Analogical analysis of the mechanical system. MM, is the mass of the magnet, M, the 
mass of the copper cylinder, Mr the radiation mass, C, the compliance of the mounting 
springs, C, the compliance of cylinder and clamping rod, C; the compliance of the coupling 
system, Cr the radiation compliance, Gr the radiation leakance, Gr’ = Grw?Cy?/(Gr?+ w?Cy?), 
Cr’ =CrG;2/(Gy2+ w2Cr”) and G is the electrical system regarded as a source of energy. 


represents the essential mechanical components of the vibrator and their physical 
arrangement, as may be seen by comparison with Figure 1. In Figure 2 (b) the 
mechanical elements—masses, compliances and leakances—are abstracted, and 
the source of motional energy is represented as a generator G. ‘The mechanical 
admittance diagram is now easily constructed (Figure 2(c)), on the principle 
that all masses are considered to move relative to a certain frame of reference, 
which may be taken as the ground itself at a distance sufficiently great to be 
completely unaffected by the radiated energy. At this stage the electromechanical 
analogy emerges, and the mechanical elements are represented by their electrical 
analogues. 

It will be noted that the masses of the mounting plate and coupling bolts 
have been omitted as negligible. The mass of the copper cylinder M, has been 
included to clarify the analysis, although both this mass and the effective mass of 
the rock M, are eventually negligible (Figure 2(d)). That the compliance C, 
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of the coupling system was significant, and that, on the other hand, all leakance 
due to friction could be neglected, were deduced from the measurements to be 
described later. For convenience in the calculations the radiation leakance 
and radiation compliance have been represented in Figure 2(d) as a series 
combination, consisting of a leakance G,’ and a compliance C,’, the formulae 
for these two quantities being given in the captions to Figure 2. 


§5. THE ELECTROMECHANICAL CIRCUIT 
The complete vibrator, except for the D.c. circuits and the driving motor, is 
represented by a single electrical network in Figure 3. In its electrical aspects 
the Fessenden oscillator has been analysed as a transformer of negligible leakage 
inductance. — 
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Figure 3. Complete equivalent electrical circuit. Zm—=motional impedance, Rp=a.c. coil 
resistance, Rs=copper cylinder resistance, Rn=a.c. magnetizing circuit resistance, 
In=A.C. magnetizing circuit inductance, t=number of turns in a.c. coil, J=constant 
current supplied by the alternator and a*=electromechanical transformation factor. 


The electromechanical relations for the copper cylinder are 
Pe JOB Ld 9) fae cae (3) and U=I0SB Lug VWsene (4) 


where F' dynes is the mechanical force on the cylinder, B gauss is the magnetic 
flux-density in the air-gap, Lcm. is the circumference of the cylinder, J amp. is 
the current flowing in the cylinder, U volt is the back E.M.F., and vcm/sec. is the 
velocity of the cylinder. The equivalent electrical impedance, or motional 
impedance, of the mechanical system is thus given by 


4y,= Uj) =P Bo r=ay, sis cegeeon 


where y 1s the mechanical admittance and a2 =10-°B?L?. Hence, by incorporating 
the factor a’, in the appropriate way, in each element of the mechanical 
admittance diagram of Figure 2(d), the mechanical system has been represented 
in its equivalent electrical form in Figure 3 (0). 

Numerical values have been inserted in Figure 3 (d) for the various electrical 
properties of the vibrator that are substantially independent of frequency. 
Laboratory measurements showed that the air-gap flux-density approached a 
saturation value of 7,250 gauss for an energizing current of 3-Oamp. D.c. Since 
the circumference of the copper cylinder is 64cm., the value of the electro- 
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mechanical transformation factor a? is found to be 215c.c.s. units. By 
open-circuited transformer tests the parallel branch of Figure 3 (c) representing 
the A.C. magnetizing circuit was found to have an impedance of the order of 
100(1+j)ohms; this branch may therefore be regarded as an open circuit, and 
has been omitted in Figure 3(d). The a.c. coil consisted of 520 turns, with a 
resistance 21-5 ohms; the resistance of the copper cylinder may be calculated as 
25-5 x 10-® ohm. 

The effects of the mechanical elements of Figure 3(c) depend essentially on 
frequency, and, as regards the radiation elements, on the properties of the rock 
as well. From Figure 3(c) two general features of the mechanical system may, 
however, be appreciated. First, since the maximum current is desired through 
the radiation elements, the compliances C, and C, should be as small as possible— 
that is, the copper cylinder should be rigid with the rock. By using the stoutest 
possible members to clamp the cylinder, as described above, the compliance C, 
was reduced to an unavoidable minimum, comprising the compliance of the 
clamping rod together with that of the cylinder itself. From the dimensions of 
these members, and Young’s moduli for steel and copper, the compliances may 
be calculated as 24-0x10-"cm/dyne and 7-0 x10-%cm/dyne, respectively ; 
altogether, C, was thus about 31:0x10-%cm/dyne. The compliance C, 
of the coupling system will be shown later to have had the very low value 
0-9 x 10-% cm/dyne. 

Secondly, Figure 3(c) shows that the compliance C, of the mounting springs 
appears in parallel with the branch containing the radiation elements, and 
therefore largely determines the performance of the vibrator. ‘There are three 
cases of interest. At one extreme, as is obvious mechanically, if these springs 
were very stiff they would ‘short circuit’ the rock altogether. At the other 
extreme very compliant springs would admit practically no force, and the force 
applied to the ground would be that exerted by the cylinder. In between these 
extremes there is a range of values of C, for which the reactances of the two 
branches are of comparable size. For such values force will be applied to the 
ground through the mounting springs, and advantage may be taken of mechanical 
resonance to increase the force. 

Two sets of mounting springs were used in the experiments. One set, 
consisting of three coil springs of parallel compliance 2:4 x 10-®cm/dyne, 
provided a virtual open circuit in Figure 3(c). The other set consisted of three 
steel rods, each 3in. long and #in. in diameter; the parallel compliance of these 
rods may be calculated as 4-4 x 10-¥cm/dyne, giving a reactance comparable 
with that of the branch containing the radiation elements. With the exception 
of the radiation elements, which in most conditions were relatively small, the 
mechanical properties of the vibrator are evaluated, in terms of the equivalent 
electrical properties, in Figure 3 (d). 

A third general feature of the vibrator is now apparent, that, largely as a 
result of the rigid construction, the equivalent electrical impedance of the 
mechanical side is of a smaller order of magnitude than the impedance of the 
electrical side. At the frequency 600c/s., for example, the reactance of the 
600-Ib. magnet mass M, is —0-06johm, whereas the resistance R, of the a.c. 
coil is 21:5ohms. This contrast gives some measure of the difficulty of 
communicating the motion of a vibrator to rock and, in particular, of the very 
low efficiency to be expected of such a machine. For purposes of analysis it is 
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concluded that the electrical side will appear from the mechanical side, in Figure 3, 
as a source of constant current. 


§6. THE VIBRATOR IN CHALK 

The experiments to be described were carriéd out in an old chalk mine at 
Chislehurst, Kent, popularly known as the Chislehurst Caves. From 
observations of direct waves the velocities of propagation in the chalk were 
found to be V,=7-0ft/msec. and V,=4-1ft/msec. The density of the chalk 
is about 1:8gm/cm?. The effective radius of the vibrator was taken as 20cm., 
the actual radius to the centres of the coupling bolts. From these values the 
radiation impedance given by equation (1) has been calculated for the frequencies 
300c/s. and 600c/s. Table 1 gives the values of the radiation elements defined 
by equation (2) and of the corresponding series elements of Figure 2 (d). 


Table 1. Values of Radiation Elements in the Chalk at Chislehurst (C.G.s. units) 


f (c/s.) Ge wM, 1/wC, Ge eke 
300 5-2x 10? 2-8 x 10? 12-9 x 10° 3:18x10- 0-41 x 1018 
600 20-8 x 107 5-7x 107 6-4 x 10° 508x108 0-41 x 10-18 


Here the positive reactive term w M, is clearly negligible. 

The equivalent electrical properties of the vibrator at Chislehurst are fully 
evaluated, for the frequencies 300c/s. and 600c/s., in Figure 4(a) and (6), 
respectively. It will be noticed that the equivalent radiation resistance increases 
rapidly with frequency as the dimensions of the vibrator become more 
comparable with the wavelength of the radiated energy. 


P,,— 45 arms. —> P 7:3 a.nm.s.> 


Relay 
Switch 
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(a) 300¢/S. circuit. (b) 600 c/s. circuit. 


Figure 4. Equivalent electrical circuit in chalk. 


Figure 4 shows that, provided the coupling compliance is kept small enough, 
mechanical resonance may be achieved by a proper choice of mounting 
compliance. ‘T’o obtain a large magnification, however, it would be necessary to 
reduce the magnet mass; or a similar effect could be gained by mounting the 
copper cylinder on springs and fixing the magnet mass to the rock. But such a 
modification would increase the selectivity of the apparatus at the expense of 
its ability to produce a well-formed pulse of vibrations. Further, in proportion 
as the moving mass is diminished, its amplitude of vibration for a given force 
is increased, and with it the intensity of sound radiated into the air. Even with 
the minute displacements occurring in the present apparatus, special care was 
sometimes necessary, particularly in working with continuous vibrations, to 
avoid airborne interference in the receiver. 

The alternator, which was rated to deliver 2kw. at the frequency 500¢c/s., 
with an output voltage of 180 v.r.m.s., was driven by a 3-H.P. motor. Frequencies 
from 300 c/s. to 1,000c/s. were obtained by varying the size of pulley on the 
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alternator shaft; at each frequency suitable condensers were inserted in the 
A.C. circuit, as shown in Figure 4. 

The a.c. circuit could be closed, for precisely 20 msec., by a relay switch 
actuated from the receiver. The pulse of voltage thus applied to the vibrator, 
measured for convenience between the points marked P in Figure 4, is shown, 
for the frequencies 300c/s. and 1,000 c/s., in the oscillograph records of Figure 5 
(Plate*). In the photographs the lower time-base carries an oscillation of 
frequency 50c/s. as a time-scale. In contrast with the form of an explosive 
impulse, the pulse is seen to consist of substantially sinusoidal oscillations of 
equal amplitude. (The high-frequency distortion produced by the alternator is 
of no importance.) At the same time, both the onset and the end of the pulse 
are abrupt. By measurements from these records the current flowing during the 
pulse has beencalculated and is given in Figure 4. 


$7. . THE RECEIVED VIBRATIONS 

Vibrations in the rock were detected by a Rothermel pickup, Model VP-5, 
which is essentially an inertia-operated rochelle-salt crystal microphone. 
Constant sensitivity was obtained by using a patent clamping device, the 
“ Rawlplug ’, to make the microphone fast to the rock. The microphone voltage 
was amplified about 10* times by a two-stage circuit of low noise characteristic, 
and a further 900 times by the oscillograph amplifier. The oscillograph was a 
Cossor Double-Beam Model 339, with camera attachment. Comparing the 
oscillograph deflection with the ground displacement, the maximum available 
magnification in the receiver was about 10°. In the records the amplitude of 
deflection was adjusted to a suitable level by means of a calibrated gain-control. 

The time-base of the oscillograph was non-repetitive, and was triggered by 
operating a switch which at the same time actuated the relay switch in the 
vibrator circuit. Thus the initiation of the radiated pulse and commencement of 
the time-base were connected in a constant time relation; it can be shown that 
the onset of vibrations at the source occurred at a time 6+1msec. along the 
time-base. 

Vibrations were recorded at intervals of distance out to the limit of the mine, 
450 ft. from the vibrator. At each point of observation, and for each frequency, 
the disturbance produced by a single transmitted pulse was found to be 
protracted and complex, resulting from the arrival of a succession of 
compressional and shear vibrations, some propagated directly, others reflected 
from depths down to several hundred feet in the chalk, Some records are shown 
in Figure 6 (Plate). Figure 6(a) and (4), taken at the same point in identical 
conditions, are included to demonstrate the high degree of repeatability achieved 
in the experiments. It will be seen that the vibrations recorded on these two 
occasions were similar down to minute details, the most noticeable difference 
between the records being in the phase of the time-scale, which was random. 

The observed disturbances are analysed in terms of the travel-time of each 
separate arrival, and various geological conclusions are drawn (Evison, to be 
published elsewhere). On the basis of the travel-time analysis deductions are 
also made from an examination of the amplitude of the various arrivals. No 
detailed study of frequency effects has yet been attempted, although it is clear 
from the records that the apparatus offers access to such effects. 


* For all Plates see end of issue. 
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§8. INFERRED PROPERTIES OF THE VIBRATOR 

The finer mechanical details of the vibrator may be deduced from comparative 
measurements of the vibrations appearing at a convenient point, with the 
vibrator in various conditions of adjustment. Figure 7 shows records made 
130 ft. from the vibrator, at the frequencies 300 c/s. and 600c/s., with the magnet 
mass mounted first on the steel rods and then on the coil springs. Since only 
the amplitude of the radiated energy is to be observed, no prior analysis of the 
received vibrations into separate arrivals is necessary. For each record the 
peak-to-peak amplitude of the maximum deflection was measured, from the 
photographic negative, to the nearest quarter-millimetre. The gain of the 
amplifier was the same for the two 300c/s. records; at 600 c/s., the gain was 50% 
higher for the coil springs than for the rods. In Table 2 the relative ground 
displacement is derived for each frequency, taking account of the relative 
receiver gain. 


Table 2. Performance of Vibrator with Spring Mountings, relative to 
Performance with Rod Mountings 


Frequency Relative oscilloscope Relative receiver Relative ground 
(c/s.) deflection voltage gain displacement 
300 16°75 mm/15-5 mm.=1-08 1-00 1-08 
600 17:5 mm/16°5 mm.=1-06 1-50 0-71 


According to Figure 4, the ground displacement at the vibrator is proportional 
to the current flowing in the substantially capacitive branch of a parallel circuit. 
As a first approximation it will be assumed that this branch contains only the 
magnet mass M,, for the radiation elements are small, and the coupling 
compliance C, was revealed only by the present calculations. For simplicity 
the transformation factor a? will be omitted from the calculations, which may 
then be referred to Figure 2(d), the generator G being regarded as a source of 
constant current. 

If Z be the equivalent parallel impedance of the rod compliance C, and the 
mass M,, and J the current in the arm containing M, and the radiation elements, 
and if the corresponding quantities for the circuit involving the spring compliance 
be 25°C, anda] ther 


['/L=2Z'/Z=(1fjoCy+joM,)/(lfjwoCy’+joM,). — ...... (6) 
But C,’ is known to be very large, so that, simplifying, 
F]f=1 Tet CMe Se eee (7) 


It may be assumed that the ratio J’/J is given numerically by the ratio of 
ground displacements in Table 2; but the sign of the latter ratio is left 
undetermined by the measurements. Rewriting equation (7) in the form 

1/C,M,=07(1-I'/D), Peele 
it is easily shown that, for the measurements to be approximately consistent, 
I’/I is to be taken as negative at 300c/s. and positive at 600c/s. The two values of 
1/C,M, are then 7-4x10® and 41x10 respectively. Evidently the minor 
mechanical elements, so far neglected, are not without significance to the 
measurements of ‘Table 2. 

The possibility of significant resistive elements appearing in the circuit, to 
represent frictional or mechanical hysteresis losses in the vibrator, or a very 
much greater radiation of energy, has been investigated. The appropriate 


An Electromechanical Source of Elastic Waves in the Ground 321 


formulae developed for these cases have led to impossible values for the 
various elements if they are to satisfy the above observations as well. It follows 
that the dissipation of energy in the form of heat has negligible effect upon the 
mechanical operation of the vibrator. 

The simplest circuit that can be postulated, consistent with the measurements, 
includes a compliance Cy, in series with the mass M,. For such a circuit, 


SNCS Cia ey ahae a 1 


Inserting the known values of J’/J, w and M,, two simultaneous equations are 
obtained : | 


2:08C, + Cu => 10-4 x 10;3°, at 300 c/s. 
and 0-29C,+Cy= 2:6 x 10-43, at 600c/s. 


Hence Cj, the compliance of the steel rods, is found to be 4-4 x 10-18 cm/dyne, the 
same value as that already calculated from the Young’s modulus and dimensions 
of the rods; it is thus clear that the only significant compliance in the mounting 
arrangement is that of the rods themselves. Completing the solution, Cy is found 
to be 1-3x10-%cm/dyne. Reference to Figure 2 shows that this compliance 
represents the sum of the effective compliance of the rock C,’ and the 
compliance C; of the coupling system. The theoretical value of C,’ is given 
in Table 1 as 0-41 x 10-cm/dyne; the corresponding value of C3, namely 
0-9 x 10° cm/dyne, will be largely accounted for by the coupling bolts. 

The calculations show that a satisfactory degree of rigidity was achieved in 
the coupling system; there was evidently no slip in the metal joints between the 
coupling bolts and the concrete or between the concrete and the rock. Upper 
limits have also been established for the values of the radiation elements, and to 
this extent formula (1) is verified. Further, particularly in giving an independent 
derivation of the compliance of the mounting rods, the calculations have shown 
the mechanical analysis of the vibrator to be sound. It is of wider importance that 
precise information has been gained from measurements of the amplitude of 
radiated vibrations, for such measurements should also be capabie of providing 
data on propagational phenomena, and on the geological formations traversed, 
with similar accuracy. 

§9. POWER AND EFFICIENCY 

From the relative recorded amplitudes of the direct waves and certain 
reflected waves it was found that the radiation pattern of the vibrator in chalk 
was consistent with that shown in Figure 8, for which there is some theoretical 
basis. he vertical pattern for the shear energy is a pair of equal semicircles 
with diameters at the level of the floor of the mine, and, for the compressional 


S 


———__ Compressional mode 
—-—-. Shear mode 


Figure 8. Vertical radiation pattern of the vibrator, showing loci of equal displacement amplitudes. 
The continuous curve refers to the compressional mode and the broken curve to the shear mode. 
S is a point on the surface at the centre of the vibrator. 
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energy, a circle of radius one-half that of the semicircles, touching the floor of 


the mine at the vibrator. The three-dimensional pattern is obtained by revolving 


Figure 8 about the vertical axis of symmetry. Integrating the absolute amplitude, 
as measured, over the radiation pattern, the power radiated during the pulse may 
be calculated. At 600c/s., with rod-mountings, the power is found to be 
0-012 watt in the compressional mode and 0-060 watt in the shear mode, making 


0-072 watt in all. The proportion of energy radiated in the compressional mode 


is thus about 17 percent. As would be expected, this value is rather greater than 
the theoretical value for an oscillating sphere, which is about 11 per cent for the 
condition 2z7f/V.,=0-35, as in the present case (Bradfield 1947). 

The total radiated power may be independently calculated from the theoretical 
equivalent radiation resistance and the measured current supplied from the 
alternator, as shown in Figure 4(4). Taking account of the current magnification 
effected in the mechanical circuit, the power so calculated is 0-031 watt. The 
discrepancy between the two values is less than 4 db., lending support both to the 
measurements of amplitude and to the use of equation (1) as an expression of 
the radiation impedance. 

The total power dissipated in the apparatus during the 600c/s. pulse was 
about 2kw. Thus the efficiency of the vibrator as a radiator has the extremely 
low value of about 0-025 x 10-°._ The efficiency of a dynamite explosion near the 
surface of the ground is of the order 1-0 x 10-° (Willmore 1949). It is worth 
noting, however, that the energy radiated from the vibrator is in precisely the 
form desired, whereas only a small fraction of the elastic energy generated by an 
explosive source is turned to use. 


$10. CONCLUSION 

The vibrator radiates elastic energy in rock in a form that can be accurately 
controlled as to frequericy, amplitude and duration. The power radiated is small, 
_.but sufficient to be detected, in favourable conditions, through many hundreds 
of feet of rock. Much greater power could be obtained without difficulty, by 
increasing the area of coupling to the rock, by reducing the electrical resistance 
of the vibrator coil and by using an alternator of higher voltage. In providing a 
more versatile means of studying the phenomena of elastic propagation in rocks 
the apparatus should assist fundamental research, particularly within the shallower 
strata of the crust, the domain of mining and engineering geology. 
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ABSTRACT. The suitability of square plates, vibrating in contour modes, for the 
determination of piezoelectric coefficients is considered. Piezoelectric measurements on 
such plates provide checks for the theoretical solution of the various modes and in particular 
for the distribution of displacement. The recently published solutions for square plates 
form a basis for this work and the reliability of these solutions is discussed. The materials 
used were sodium chlorate, sodium bromate, quartz, ammonium dihydrogen phosphate and 
ethylene diamine tartrate (EDT). New determinations of the elastic and piezoelectric 
coefficients of these materials have been made, with the exception of EDT, the data for 
which were published in a previous paper. 


§1. INTRODUCTION 

URING recent years, the vibrations of crystalline bodies have been 
A) widely studied. Rough approximations for the solution of free vibrating 
plates of anisotropic materials were obtained by Mason (1934) for the 
face shear mode of square plates and by the author (Bechmann 1941, 1942) for 
the contour modes of vibration of square and rectangular plates. More 
systematic studies made by Ekstein (1944), Mahly (1945), Mahly and Trésch (1947) 
on thin square plates have resulted in better approximations for the solution of 
both the frequency and the distribution function of the displacement. The 
author’s (1941) approximate solution describes the shear mode and two 
longitudinal modes of vibration; a third observed longitudinal mode (Bechmann 
1942) (the frequency constant denoted by N,) was not explained. The solution 

of this mode has been derived by Ekstein (1944). 

The form of vibration can be studied from lycopodium patterns and by various 
optical methods. The solutions can also be verified from results of electrical 
measurements of frequency and of the elements of the equivalent electric circuit 
of plates with known elastic and piezoelectric constants, since the elements of 
the equivalent electric circuit are functions of the distribution of displacement. 
Once the solutions have been so checked, they can be used to derive, dynamically, 
unknown constants of other materials. In some cases, the same coefficient of a 
crystal can be obtained from different modes. For example, the piezoelectric 
coefficient d,, of a crystal belonging to a cubic class can be determined from four 
different modes of square plates and also from the longitudinal mode of narrow 
bars, thus giving five independent results and forming a check of the theory on 
which the formulae were based. Differences in these values, after allowing for 
errors in measurements, are due to the approximate nature of the solution. 

With the recent extended use of synthetic crystals as piezoelectric resonators, 
one of the problems associated with the investigation is the production of cuts 


with small temperature coefficients of frequency and of the elements of the 
This necessitates a knowledge of elastic and 


equivalent electric circuit. 
Ne 
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piezoelectric constants and their behaviour with temperature. Some of these 
constants can be determined using the results on the various modes of vibrating 
square plates mentioned above. 


§2. CONTOUR MODES OF VIBRATION OF THIN SQUARE PLATES 
AND LONGITUDINAL MODE OF NARROW BARS 


In the following, using a Cartesian coordinate system, a thin square plate of 
the thickness ¢ is assumed to have its normal parallel to the Z axis and its sides 
of length / parallel to the X and Y axes, with the centre of the plate at the origin. 
The axes of the coordinate system do not necessarily coincide with any of the 
natural axes of a crystal but the usual ‘ dashes’ denoting arbitrary orientations 
are omitted. The three longitudinal modes of vibration will be denoted 
modes 1, 2, 3 as introduced by Ekstein (1944), who describes them in detail with 
diagrams. For comparison, the equations for a narrow bar are included, where 
a narrow bar is defined as having its length / along the X axis, its breadth 6 along 
the Y axis and its thickness ¢ along the Z axis, the centre at the origin of the 
coordinate system. 

The solutions for the different modes of vibration for free and undamped 
plates vibrating in the fundamental frequencies are considered. The conditions 
of these plates are X,=0, Y,=0, Z,=0 on the surfaces, and also throughout the 
plate (Voigt 1910, Ekstein 1944) and the strain—-stress relations reduce to the 
equations 


—Xy =SyX_+5y2Vy+ eX y; 
=p = SaX BE Sead, Hoag, po MODOG> Tae (1) 
= Hy =SygX qt Seg Vy + SegX y- | 
X, Y,, X, are the average values of the stresses taken over the thickness of the 
plate. Solving equations (1) by determinants: 
Xp =VuXet+ VieVy t+ VreXy> 
Meo are CT te at ee (1a) 
Xy =YVigX2 t+ Yo6eVy + YeoXy» 


where y;, =x; are elastic moduli (stiffness factors) first derived by Voigt (1910). 
For*quasi-isotropic conditionsin,; =s74)3)595 0) Sig 0ieey eels Seen (2) 
the expressions for y,,, become: 


S44 E 
Ye 99 ye ees 
C Spy Syghiy lil i0r 
Sie Eo (3) 
112 i 2 3 > Soler 6h POMveNre ers 
Sheers Te 
1 
Yee=—> Yie=Y2e=9, 
566 
where E=1/s,, is Young’s modulus and o = — 545/8;; =y4/y1, is Poisson’s ratio. 


For the more general conditions, when 


S11 =Sgq AN S46 = So6 
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these expressions are: 
Hy = Wyap = 544556 — $16" 
Iyy2 = S46” — $12566> 
Pe si ee ee ee (5) 
Hy16 = Wyeg =(s12— S11) $165 


TL = (541° — 497)5¢¢— 2(811 — Sia) S16? 

Under conditions (2) and (4) the relationship (s,;—542)(y1;—712) =1 holds. 
The y,;, (7, k =1, 2, 6) are only applicable for the Z plane; for the other coordinate 
planes similar expressions hold, but the indices 7, k are changed in a cyclic order. 
For each coordinate plane, the six expressions y,; transform in a similar manner 
to the corresponding elastic moduli ¢;,. The transformation can also be made 
using the transformed s,, in the expressions (3) or (5) defining y;,.. In the more 
general case of oblique planes the second method is the simplest. 

The frequency vy =w/2z of the different modes of vibration considered can be 
calculated from the general equation, 


Perkin, (Ct) oN anindle: (6) 
where p is the density, & the wavelength constant and y the elastic modulus, which 
assumes different forms for the different modes of vibration. ‘Table 1 contains 
expressions for y for the three longitudinal modes of square plates, the longitudinal 
mode of narrow bars and the contour shear mode of square plates, all under 
conditions (2). All values for the elastic coefficients s,;, used to calculate the 
expressions for y are the isagric values sf. The expressions for the components 
of displacement uw, v for these modes are also included.* For the longitudinal 
modes k=7//, and for the shear mode k is defined by equation (8). A;(z=1...5) 
is a small arbitrary constant expressing the amplitude of displacement. 

The following remarks may be made with respect to the different modes: 

(i) Longitudinal mode 1 of square plates. Identical solutions have been derived 
by different methods by Petrzilka (1935), Bechmann (1941) and Ekstein (1944). 
The solution is exact for plates fulfilling conditions (2) and the more general 
conditions (4). 

(ii) Longitudinal mode 2 of square plates. 'The solution for this mode was 
first derived by Ekstein (1944) for conditions (2), but is valid for both 
conditions (2) and (4). The expression for y corresponds with that given by 
Petrzilka (1935) for his modes A and B, but his equations for the displacement 
do not agree with Ekstein’s. Increasing differences between observed values 
for frequency and values calculated with y given in Table 1 are found for increasing 
values of T=(749?+7162)"2/y1. ($5.2 (ii)).. In the case of cuts fulfilling 
conditions (2), tT =o. . 

(iii) Longitudinal mode 3 of square plates. 'Vhe solution quoted in Table 1 
was given by Ekstein (1944). An isolated mode can only be obtained under 
conditions (2) or under the more general conditions when $1; =5g2 and S25 = — 546, 
and there is no coupling with the shear mode. In the author’s earlier 
approximation (Bechmann 1942), the frequency is determined by the elastic 
modulus 


Sh OCR ee ee ree (7) 


* The argument e’” is omitted from Tables 1 and 2. 
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Table 1. Solution of various Modes of Vibration 


y Components of displacement 

Longitudinal = u=A, sin kx cos ky 

mode 1 Yi Y12 v=—A, cos kx sin ky 
Longitudinal u=A,(cos ky — 47) sin kx 

mode 2 Yat v= —A,(cos kx — 47) sin ky 
Longitudinal 8 u=A, sin kx 

mode 3 Yut panies v=A, sin ky 
Longitudinal 1 

mode — u=A, sin kx 

narrow bar a4 

Contour u= A, sin RY ..< 
shear mode Ys v=A;sinkx+... 


The more rigorous expression given in Table 1 can be written as: 


race ea hk 


This differs from the expression (7) only by the term in square brackets which 
depends on a, but the author’s formula for the displacement has a different form 
from Ekstein’s. Mahly (1945) criticized Ekstein’s formula for y and also 
formula (7). He mentioned a better approximation for the solution of this mode, 
which has not been published, and gave his results for the frequency in the form 
of a curve plotted against o. Also on the same graph are curves of results derived 
from Ekstein’s expression for y (Table 1) and from (7) providing a comparison 
of the accuracy of these solutions. In addition, he gave a graph showing the 
distribution of the displacement for o=0-3, which differs slightly from the 
distribution given in Table 1. 

(iv) Longitudinal mode of narrow bars. Line 4 of Table 1 contains the 
corresponding expressions for longitudinally vibrating narrow bars, which are 
included for comparison and for the determination of the appropriate elastic 
coefficients 5;,. 

(v) Contour shear mode of square plates. Approximate solutions for shear 
vibrating square plates have been given by Mason (1934), Bechmann (1941), 
and Ekstein (1944). Recently a closer approximation has been derived by Mahly 
and 'T'résch (1947) which takes into account longitudinal effects in order to fulfil 
the boundary conditions. The expression for y and the main terms of the 
expressions for the components of displacement * are given in Table 1, but the 
terms covering the longitudinal effects are omitted as they only influence areas 
near the edges and do not affect the considerations in the following sections. 
The wavelength constant for this mode is 


k=2k/l Aces (8) 
where tan et K=0.05 fae oe ere (9) 


The first root of equation (9) isx=2-0288. But to determine the frequency from 
equation (6), a corrected expression for k should be used: 
k=2kA/l weeenn( 10) 
1 /x?-2\8? s 
h A=1- —(—.— ue =, patio 
where 1 Fe (5 i. 5) pes ae ra (10 a) 
Mahly and Trosch discuss the accuracy of the frequency solution in their paper. 


“ The author has been informed of the displacement expression by a private communication 
from Mahly, and details will be seen in his Thesis (Ziirich). 
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§3. ELEMENTS OF THE EQUIVALENT ELECTRIC CIRCUIT 

The properties of any piezoelectrically excited mode of a lightly damped 
mechanically vibrating system having many degrees of freedom can be represented, 
near resonance of the operating frequency, by an electric circuit consisting of 
a series combination of capacitance, inductance, and resistance, shunted by a 
capacitance approximately equal to the static capacitance of the electrode 
arrangement. ‘The dynamic capacitance and inductance are the stiffness and 
effective mass of the mechanical oscillator transformed into electrical terms by 
the piezoelectric effect. Since the mechanical and electrical systems are 
equivalent, the energy at any instant due to the piezoelectric excitation is equal 
to that of the mechanical system. This relationship is given by equations for 
frequency, and kinetic, potential and dissipated energy, but only those for 
frequency and kinetic energy T are used: 


wee =k (2) = som ae oct Gall 1) 
=2rv= ; PAG By at iea od armoire 0 ( 


f 3:feu\? fea\* 1 
ak a 


the integral being taken over the volume V of the vibrating plate, where L,, andC,, 
_are the inductance and capacitance of the equivalent electric circuit and J,, the 
total current due to the piezoelectric polarization. 
From equations (11) and (12), the inductance and capacitance are 


ar I 


It is assumed that the faces of the plates normal to the Z axis are covered with 
a thin conducting layer. ‘The component of the piezoelectric polarization P, 
is produced by the mechanical stress and, therefore, the current J, is obtained 
from the integral over the covered area a, 
eg $4 [ 7 (fla Ve dgX \aa oe (15) 
Primes are introduced here for the piezoelectric coefficients to emphasize the 
use of an arbitrary coordinate system. Equation (15) can also be used to study 
the effects of partial plating. 

Since piezoelectric crystals of practical interest are generally lightly damped; 
the solutions for free and undamped vibrations, given in § 2, can be used without 
introducing any appreciable error. 

The following assumptions for the piezoelectric coefhcients are made for the 
various modes considered in §2: 


longitudinal modes 1 and 2: d=d3;'= — dy’, 34 =0 

longitudinal mode 3: d = dsy' = dso’, d3_ =9 eae (16) 
longitudinal mode for bars: d=d;,' 0, dg’ =0, dz =0 

shear mode: d=dz, £0. ds;'=0, dso =0 


Y 


All other piezoelectric coefficients are assumed zero. The primed symbols with 
the same suffixes for the piezoelectric coefficient are different for the various modes. 
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Table 2. Energy, Current and the Constant 8 for various Modes of Vibration 


fh Ip B 
Longitudinal 4? w 512 
es 1 re pw*V 8Ajyd— Ee 
Longitudinal Ae rel. or cis 512 a ne 
mode 2 Aovern rs ce bie’ 8Asyd> 1 8 m\8 
Longitudinal A;?* w 64 
cago! Aaya te 2 
Longitudinal Az . 32 
mode Se pw?V 2Ayydwb = 
narrow bar a 
Cont A;* 3 ’ 47? 
ey 5 pw?(1+cos*«)V 4A.yduwl sin k RD) 


From equations (12) and (15), using the formulae given in Table 1, the 
expressions for kinetic energy T and electric current J, can be calculated, since 
X,=yx,;. The results, assuming the faces are fully plated, are given in Table 2, 
where y and d are abbreviations for the resultant elastic moduli and piezoelectric 
coefficients as defined in Table 1 and in equations (16). The dynamic 
capacitance C, can be obtained from J and J,, using equation (14). By 
eliminating y, using equation (6), and introducing the capacitance constant 
I =C,,t/a, and the frequency constant N=vi, a general expression is obtained 
for the capacitance constant, which holds for all the cases considered: 


T= pBd*N?,: fy > Oe (17) 
where 6 is a constant, dependent only on the distribution of displacement, having 
the form given in Table 2 for fully plated crystals.* Equation (17) also holds 


for a first approximation of the effects of partial plating and £, then, depends 
on the arrangement of the electrodes. Solving equation (17) for d: 


fe (| m) aan (18) 


‘The frequency constant N is obtained from equation (6), 


1/2 
ie ”)  stanitesnice OEE (19) 


where y 1s taken from ‘Table 1 and & is as previously defined. 


For the shear mode, if the earlier solution for the distribution of displacement 
given by Mason (1934) and the author (1941): 


u=coskxsinky, v=sinkxcosky (k=7/l) 


is used, then £,=128/7* and, comparing with the value in Table 2, the ratio 
B/B, =7°/32«>(«? +2) =1-19; substituting the value of 8, in equation (18), values 
of d about 9-27%, higher would be obtained. Similarly, if the author’s (1941) 
distribution formula for mode 3: w=sinkxcosky, v=coskxsinky is used, 
then 8,=512/7* and, comparing with Table 2, the ratio 8/8, =7?/8 =1-23, 
which results in a value for d which is 11-1°% higher. 

_* Equation (17) can also be applied to the case of flexurally vibrating narrow bars with two 
pairs of electrodes. For a bar with identical orientation to that of the longitudinal bar considered 


(d=ds;’): N=vl*/b. For the fundamental mode, B=(97? tan%e)/e®, where tan e-+tanh e=0 
(«=2-8650). 
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The solution of any mode comprises the distribution of the displacement 
and the resultant elastic modulus y. If measured values of N and I are 
substituted in equation (18), the value of the piezoelectric coefficient d depends 
only on f and, provided d is known, the distribution of displacement can be 
checked. Similarly, equation (19) tests the accuracy of the approximation of y, 
if the appropriate values of s,, are known. Since equations (18) and (19) are only 
related by N, if the measured value of N is used, equation (18) still holds, even 
if the approximation of the solution for y is not very close, or if the solution is 
affected by piezoelectric secondary effects. 


§4. SELECTION OF SPECIMENS USED EXPERIMENTALLY 

In order to check experimentally the formulae for the various modes on 
different crystals, particular orientations of plates and bars must be chosen, which 
fulfil the assumed conditions. Table 3 contains the orientations of plates and 
bars of sodium chlorate, sodium bromate, quartz, ammonium dihydrogen 
phosphate (ADP) and ethylene diamine tartrate (EDT), which were measured. 
Isolated vibrations of mode 3, excited by a field parallel to the thickness of a 
plate fulfilling the conditions (2) and (16), can only be obtained with crystals 
belonging to the cubic system. 


Table 3. Materials and Specimens used for Measurement 


Crystal system Cubic Trigonal Tetragonal Monoclinic 
Symmetry symbol Ay Ds Va C, 
Piezoelectric coefficient dia dane dis 36 dis 36 
Long. mode 1, 2 plate Xo°, 45° Yo, 45° Xoo 45° Zoe, 45° Xo, 5° Zoe, 45° 
Long. mode 3 plate P4s° 55° — — = 
Long. mode bar Dy Ky xy Ly hp Zy 
Shear mode plate Xoo 0° Yo, 0° Xoe o° Zo°, 0° Xo,9° Zoe 0° 
Crystal measured Sos Baca quartz ADP EDT 


(i) Explanation of the symbols of specimens used in Table 3. Square 
plates Xp» 4» Zo: 9, refer to plates, perpendicular to the X, Z axis with sides 
parallel to the other axes; square plates Xo» gs, Zo», 45° are similar plates to the 
above with their diagonals parallel to the other axes. ‘The symbols Y4 9°, Yo, 45: 
refer to square plates similar to the above, but rotated about the X axis through 
an angle 6 in the notation of Bechmann (1942). The square plates P45, 55. (used 
in the cubic system only) are plates which lie in the crystallographic plane (111), 
having one side in the coordinate plane. Narrow bars are denoted by X,, Zy, 
meaning that the length of the bar is in a plane perpendicular to the X, Z axis 
and rotated about this axis through an angle 7, measured in an anticlockwise 
direction. ; 

(ii) Elastic and piezoelectric coefficients of the specimens. . The elastic and 
piezoelectric coefficients, determining the vibrational properties of the defined 
square plates, belonging to the cubic, tetragonal and monoclinic systems are 
given in Table 4 for shear modes and in Table 5 for longitudinal modes 1 and 2. 
In Table 5 the elastic coefficient s,,', which equals sy’, is obtained using the 
positive sign, and 5,,’ using the negative sign in the expressions in the third 
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column. Thus it is obvious that’s,,’—51.’=4s, (¢=4 or 6). Primed symbols, 
used in the headings of Tables 4 and 5, distinguish the rotated cuts from the cut 
Ze, o, On which Table 1 was based. For the rotated cuts, the symbols of the 
elastic constants in Table 1 should be read as primed and the expressions in 
Tables 4 and 5 are the corresponding transformation equations. The 
transformation equations for the plate P,;. 55. in the cubic system are: 


Sir’ = So’ = $[2(511 + Syn) + Sqa], $16, =Sy6 =9, 
Sy’ = $512 + 75 [2(511 + S42) — Sua] d= ds)’ = dye’ = — V/3 dy/6, 5 «++ (20) 
See = $ (S11 — S12 + See)> 36 = 0. 


Table 4: Transformation for Elastic and Piezoelectric Coefficients 
for Shear Mode 


Symmetry Cut Sac Soo. S42) Sea a=dee 
‘te Xo°, 0° Sis Sia Sig S44 dys 
Vv Xo, 0° Sui $33 Si3 S44 dy4 
d 7, d 
(020° Sia Si Sie S66 36 
G Xo°, 0° See S33 S23 S44 dys 
0°,0° Sui S22 Si2 S66 d36 


, , - ea Fe 
S16 =Seg =0; d31’=d3_’=0 


Table 5. Transformation for Elastic and Piezoelectric Coefficients for 
Longitudinal Modes 1 and 2 


$1, using+sign , d 


Symmetr ,Cut ; ‘ s Sie sone : 
y Ba te ’ $49’ using —sign es BS * >, dsj ade 
1 
rT Xo°, 45° $12(S11 +512) +544] 2(S11—S12) 0.,; aay 
ote L a | 1 1 
Va S Xo » 45° H[811 +533 + 2513 +54,] S11 +5S33—2813 $(S33—S11) tare 3d 14 
; : : ae 
\ Zoe, 45° #[2(S11 +S) 2) + S66] 2(S11—S12) 0 336 
Y " a | ; ‘ = 
cs Xo,° 45° t[Seo+S33 + 2Sog3tS4a] Sgat+S33—2823  3(S33—Saa) 314 
ies 1 if ! eed ; 
Oba ane ESy1 See 289+ Se] - Sy, Sa2—2819' $(S22—S11) 3436 
=) 


The long expressions for the elastic coefficients of plates Yy 9, Yo, 45> belonging 
to the trigonal system, symmetry class D, have not been included as they can be 
found in the author’s paper (1942), but the piezoelectric coefficients as function 
of the angle @ in equations (16) are 
shear mode: | 
d=d3, = —2d,, sin @ cos 0+ d,,4sin? 0, 
longitudinal modes 1 and 2: 
d= ds)'= — dgy' = —d,, sin8cos0+4d,,sin?6. | 


‘The formulae for s,,’ for longitudinally vibrating bars are well known; the values 
for s,, for bars X,;. and Z4;. for the crystal classes T, V, and C, are identical 
with those of plates Xq. 45. and Zos 45. as given in Table 5. The piezoelectric 
coefficients for these classes are 

X, cut:d=ds,'=d,,singcos; Z,,cut:d=dsy' =dsgsinyscos ......(22) 
and for class D, 


X, cut: d=d;,'= —d,, cos?s+d,, sin’ cos yz. 3: at(Z9) 
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The values for y;, can be determined from Table 5, for example, these values 


for plates Zp. 45» symmetry class V, (conditions (2)) can be written using 
€quations (3) 


Yu =Y22' = eekly a es! : | 
oS gat eee — So Yonw arene en 5a 
Z2SutSig Seg’ oe 2(541— S12) (24) 
i i i . f PIs 
ein me HE nad S00. eth eee | 
and for plates Zp>, 45. symmetry class C, (conditions (4)) using equations (5) 
Mi). eS ty 28 japon LS 2 eetes LO ego 
4 $4599 — S49” Seo’ Mea 4 $1599 — 549” See” 
p ae 1 543 +599 +2545 : in Veter 
66 = - —-C_ — — na ° 
4 544592 — S49" Eas 4 544899 — $49” 
batt. (25) 


Similar expressions can be obtained for the other cuts mentioned in Table 5. 


§5. EXPERIMENTAL RESULTS 

Measurements of frequency and capacitance were made on a large number of 
specimens, as described in §4, due heed being paid to remarks in § 6 of the author’s 
previous paper (Bechmann 1950), which deals with the experimental technique. 
All specimens were silver plated on the main faces by evaporation and were held 
by spring clamps at their nodal points. The accuracy of the measurements is 
within the order of 0-2% for the capacitance and 0-02% for the frequency. All 
values given are average values based on sets of measurements on several 
specimens of various dimensions for each cut and are adjusted for 20°c. 

The dimensions and units of the quantities in the following tables are: 
frequency constant N in kc/s.mm., elastic coefficient s,, in 10-%cm? dyne~}, 
elastic modulus y;, in 10'°dynecm?, piezoelectric coefficient d;, in 10-° c.G.s., 
E.S.U., and temperature coefficient of d,,, (Td,;,) in 10~4/°c. 


5.1. Piezoelectric Coefficients 


The piezoelectric coefficients derived from these measurements by formula (18) 
using the values for 8 in Table 2 are given in Table 6 which also includes the 
values for the density p in gm. cm”?. 


Table 6. Values of Piezoelectric Coefficients d;, in 10-8 E.s.u. dyne~* derived 
from Different Modes and their Temperature Coefficients T'd,;, in 10~4/°c. 


Longitudinal Shear 


Material dix Mode 1 Mode 2 Mode 3 grt! be pe Tdi, 

NaClO; 5.9 ; 19 5 Dz Sop -+- 80 

p=2-490 ais 5223 5°37 5°19 5-24 5D) 

NaBrOs; 7.4 , 75 7.68 96:5 

3-339 dia WET 121 TIE Wet] 6 

Quartz alee, 6:93 6:88 — 6°95 6:78 —2:15 

p=2:654 dy4 —2:18 S221 _ —2-20 —2-11 +-12°9 
ZA) E aie Sez 5°43 —' 5-44 5-20 —180 

p=1-803 dse —145 —149 _ —145 — 144 —28- 
EDT Say —30-1 —29°8 = —30°5 —29-9 —7°5 


PH 53811 dig 2155-4 » 565 - — 55-6 54-6 wv e165 DB 
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(i) Longitudinal mode 1 of square plates. The solution for mode 1 is exact 
and as previously mentioned (Bechmann 1950) it is thought that the correct 
values of the piezoelectric coefficients are most likely to correspond with the 
results of this mode. The values obtained from the other modes are compared 
with them. The piezoelectric coefficients d,, and d,, for quartz are derived from 
formula (21) using various combinations of different orientations. The coefficients 
derived from this mode are independent of the orientation of the specimens. The 
observed values of CT cut quartz plates (@=128°) and DT cut (@=38°) are 
ds,’ =2-685 x 10-8 and d,’= —3-775 x 10-8 respectively. The coefficients dj, 
and ds, for ADP have opposite signs as obtained from oblique cuts Xy 45°. 

(ii) Longitudinal mode 2 of square plates. The values of the piezoelectric 
coefficients from mode 2 are about 2 to 3% higher than those from mode 1 with 
the exception of those of quartz plates which are 0-7 and 1:7% lower for DT 
and CT cuts respectively. The value for 7 (defined in §2(ii)) for quartz CT cut 
is 0:07 and for DT cut 0-3; for all other measured plates 7 is between 0-35 and 0-63 
(see Table 8). The deviation appears to be slightly dependent on 7. Table 8 
shows that the frequency constant N varies greatly with 7 and differences between 
observed and calculated values for N are as large as 25°% for the measured 
specimens, but the distribution function appears to be only slightly affected by r. 
The values for the piezoelectric coefficients of quartz obtained from mode 2 
depend to a small extent on the angles used for solving the combinations. 

(ii1) Longitudinal mode 3 of square plates. The values for the piezoelectric 
coefficients from mode 3 agree very well with those from mode 1. However, 
since for both materials, c=0-25 (see Table 9) and according to Mahly (1945) the 
distribution of displacement appears to be slightly dependent on a, it cannot be 
deduced that this good agreement holds for other values of o. The results from 
the two measured crystals confirm Ekstein’s distribution formula. 

(iv) Shear mode of square plates. A comparison between the shear mode and 
mode 1 shows that the values of the piezoelectric coefficients from the shear mode 
are 0:5 to2% lower. This difference appears to be a function of 1/ and increases 
with increasing \/ (defined in §2(v)). For ADP plate Zo: 9:, wis 0-1 (see Table 11) 
and the difference between the shear mode and mode 1 is only 0:-6%. As for mode2, 
the values for the piezoelectric coefficients of quartz determined from this mode 
depend to a small extent on the angles of the specimens. 

(v) Longitudinal mode of narrow bars. As can be seen from equations (22) 
the same piezoelectric coefficient can be determined from bars at different angles. 
For bars X,;. of sodium chlorate and sodium bromate and bars Z,;- of ADP a 
complete agreement exists with mode 1. These are the values given in Table 6, 
and the figures quoted for the other piezoelectric coefficients from this mode 
are average values. ‘lhe value for d,, for quartz was obtained from narrow bars 
Xo: (bars parallel to the Y axis). Observations have shown that values for the 
piezoelectric coefficients obtained from the longitudinal mode of narrow bars 
are slightly higher than those from mode 1. 

For comparison, ‘Table 7 contains values for the piezoelectric coefficients 
obtained dynamically by Mason for sodium chlorate and bromate (1946 b), for 
quartz (1943) and for ammonium dihydrogen phosphate (1946 a) and the static 
values for these crystals measured by Spitzer (1938). 

(vi) Temperature coefficients of piezoelectric coefficients. ‘Table 6 contains the 
temperature coefficients Td;, of the piezoelectric coefficients d,;, obtained from 
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Table 7. Comparison of some Recently Published Values of Piezoelectric 
Coefficients d;, in 10-8 £.s.u. dyne 


Material d; Spitzer Mason New values 
NaClO, aia —5-23 6:1 : 5-23 
NaBrO, dis PAIL 8:2 eden 

diy 6-94 6-76 6°93 
Quartz 

di, — 2750 —2-1'8 

“thon 4°53 5 0) 5 27 

ADP 
dsz_ —136-7 148 —=145 


mode 1. They are average values in the temperature range 15 to 45°c. A 
summary of previous measurements of the temperature coefficient Td,, for quartz 
can be found in Cady (1946).* No published values for Td,, are known.t 

The dependence of the temperature coefficients Td on the angle @ for quartz 
plates Y, ,;- and on the angle y for narrow bars X,, is of practical interest as they 
determine the temperature coefficients of the elements of the equivalent electric 
circuit. The temperature coefficients, Td,,' for mode 1 of square plates derived 
from equation (21) and Td,,’ for longitudinal mode of narrow bars, derived from 
equation (23) are shown in Figures 1 and 2 respectively. Zero points occur 
at @=47° for plates Y, 4;- and at y=28° for bars X,. For the DT cut (6 =38°) 
Td3;' is —0-5 x 10-* and for the CT cut (6 =128°) Tad,,’ is —6 x 10-*. The values 
are also small for Td,,' for X,, cut quartz bars, near 4 =0. 


1 Nene} al a et 
0 2 4240 &«260 #80 100 120 (140 60 180 0 20 «+40 «6 #«80 


6 (degrees) W (degrees) 
Figure 1. Temperature coefficient of the piezo- Figure 2. Temperature coefficient of the piezo- 
electric coefficient of square quartz plates electric coefficient of marrow quartz bars 

Yp o2 and V% 45 as function of the angle 0. Xy as function of the angle #. 


5.2. Frequencies of the Modes considered 


The frequencies of the various modes calculated from expressions in $2, in 
particular from Table 1, are compared with the measured values. 

(i) Longitudinal mode 1 of square plates. From observed frequency constants 
of mode 1, the elastic shear coefficients $44, 5g, are determined and given in Table 8. 


* The value for Td,, derived from Clay and Karper (1937) should be read as —1+1 x 10-4 instead 
Ore 1x 105°) Cs 

+ Note added in proof. Since this paper was written two papers have been published (a) by 
A. C. Lynch (Proc. Phys. Soc. B, 1950, 63, 890) of this laboratory describing earlier work and (6) 
by R. V. Cook and P. G. Weissler (Phys. Rev., 1950, 80, 712) of the Bureau of Standards. ‘The 
values for Td,, derived from the data given in the latter seem rather high. 
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But for quartz as plates rotated at angles 0 =38° and 128° are used, the coefficients 
See = See COS? O + 544 sin? 6 — 451, sin 8 cos @ are given. Since the solution for mode 1 
is exact and the conditions (2) or (4) are fulfilled for the measured specimens, the 
values for s,, and sg, are regarded as correct and the following calculations for the 
other modes are based on these values. Measured frequency constants for 
mode 1 for square quartz plates Yj 4; as function of the angle 6 have been 
published by the author (Bechmann 1942, Figure 5). In this paper, the frequency 
constant of mode 1 is denoted by Nj. The variation of frequency constant with @ 
for mode 1 has been recalculated with improved values for s,,, (see Table 10) and 
is plotted in Figure 3. The earlier observed values for the frequency constant 
have been taken from the paper mentioned and are shown in Figure 3 by circles. 
The agreement is very good and the earlier discrepancy between calculated and 
observed values has vanished. At the angles @=91° and 153° the elastic 
coefficient 5,5’ is zero and changes its sign between these angles. 
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2400 
2200. | a = 
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Figure 3. Frequency constants N of modes 1 and Figure 4. Relative difference D between 
2 of square quartz plates Y, 4,0 as function observed and calculated frequency 


of the angle 0. constants of mode 2 of square 
plates as function of r. 


(11) Longitudinal mode 2 of square plates. The values for the frequency 
constant of mode 2 are calculated from the expressions given in Tables 1, 5 and 8 
using values for s,;, from ‘Table 10 and are compared with the observed values in 
Table 8. The relative difference D=(Neaje—Nops)/Nops is given in the last line. 
It has been found that this difference D increases rapidly with increasing 7 (see 
§2(ii)) as shown in Figure 4, where D is plotted against 7. The curve is an 
empirical approximation of the form D=0-5(72+74), which agrees with the 
observed values very well. The points marked by circles refer to values obtained 
from plates fulfilling conditions (2) and the crosses to values from plates fulfilling 
conditions (4). Earlier measurements of the frequency constant for square 
quartz plates Yy 45. made by the author (Bechmann 1942, Figure 5) for mode 2 
and denoted by N, are also included in Figure 3 and are shown by crosses. The 
frequency constant for mode 2 for these quartz plates has been calculated as a 
function of the angle @ with the given formula, using values of s,;, from Table 10 
and is also plotted in Figure 3. The agreement for this mode with the observed 
values is very satisfactory for small values of r (r<0-1). Between 6=90° and 160°, 
7 is less than 0-075. Between 10° and 60° there is a deviation which has its 
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maximum at about @=38° coinciding with the maximum of 7+=0-303 the 
deviation at this point is 4+-79%. The differences D for these quartz plates in the 
range 10° to 70° are included in Figure 4. Between @= 90° and 150° the frequency 
constant curves for mode 1 and 2 are interchanged, as 51.’ changes its sign 


(§ 5.2, (i)). 
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(iii) Longitudinal mode 3 of square plates. Table 9 contains the values for 
| Six’) Sig’ and o” for sodium chlorate and bromate for mode 3, derived from 
| formulae (20), together with the frequency constants calculated from 
(a) expressions for y given by Ekstein (Table 1) and (6) the author’s expression 
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(equation (7)), using the values of s,, from Table 10. Also included are the 
observed values of N and the relative differences D compared with both calculated 
values (a), (b). From Mahly’s Figure 1 (1945) these differences D are 0-95 and 
2:9% respectively, giving a good agreement with the values in Table 9. 


Table 9. Data for Longitudinal Mode 3 of Square Plates, 
sy in 10-%em*dyne. NV in ke/s. mm. 
‘ / , / N N, 
Material Neos S11 ve o () D, D, 
NaClO 3 2040 30°28 —7°51 0-248, 2060., 2100. 0:010 0-029 
NaBrO 3 1970 24°21 —6:07 0-250, 1992., 2031., 0-012 0-031 


(iv) Longitudinal mode of narrow bars. ‘The longitudinal mode of narrow 
bars was used for a new determination of the appropriate elastic coefficients 5;;,; 
the results are given in Table 10. The values of s,, and sg, were derived from 
mode 1 and taken from Table 8. A comparison of these values with the other 
recently published values would lengthen this paper unnecessarily. The values 
of s,,, of EDT are given by Bechmann (1950). 


Table 10. Summary of the New Values of Elastic Coefficients s,;, in 
10°" cm dyne* 

Material San S33 Si2 S13 Saa See Sia 
NaClO3; DSS) — O05 — 85°4 — — 
NaBrO,; 20-4 — —4-8 — 65°65 — = 
Quartz NACI 9-60 —1:79 —1-22 20:00) &29-13 —4-48 
ADP 18-10 435 ig) —11-8 11523 164-6 — 


(v) Contour shear mode of square plates. ‘The observed values for the 
frequency constant N for the defined square plates are shown in Table 11. From 
these values s,,/A? or sgg/A? can be calculated using equations (10) and (19), 
where yp is the correction (10a) stated by Mahly and Trosch. As the value of u 
can be deduced from Table 10, A can be calculated and hence s4, or sgg. These 
values agree very well with the corresponding values obtained from mode 1 given 


in Table 8, showing that Mahly’s formula gives a very close approximation to 
the solution. 


Table 11. Data for Shear Mode of Square Plates, s,; in 10-18 cm? dyne~? 
N in ke/s. mm. 
Material ; : AA? Sit 
ateria Cut Shi IN [arse sy/A be A eneirena 
NaClO, Xo°,0° Sa4 1362 90:29 0:2681 0:9740 85:66 
NaBrO3 Xo°0 S44 1340 69°56 031079 9110-97209) 865372) 
Quartz { Y 39°,0° Stee 2073 36°56 0:2984 0:9726 34-59 
128°,0 See. 3080 16°56 0°8155 0:9547 15-09 
ADP 4 Xo°,0° Saa 1383 120-9 0:2671 0:9741 114-7 
Z°,0° S66 1166 170-1 0-100 0:9841 164-8 
EDT JS Xor.0° S44 1150 205-0 0:3564 00-9701 192-9 
Zo°.0 566 1166 19OeS 0:1826 0-9786 


191-0 
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An Absolute Measurement of the Susceptibility of Tantalum 
and other Metals 


By F. E. HOARE anp J. C. WALLING 
Physics Laboratories, University of Leeds 


MS. received 25th August 1950, and in amended form 6th December 1950 


ABSTRACT. 'The magnetic susceptibility of tantalum has been determined absolutely, 
the value obtained for the mass susceptibility being yr, x 10°=0-8490+0-0006 at 20° c. 
Small specimens of tantalum have been used to determine, by a comparative method, the 
susceptibilities of specimens of platinum, palladium and rhodium at 20°c. The values obtained 
are: xpt_ Xx 10®=0-9712+0-0007, ypq x 10®°=5-231+0-004, xrn x 10°=0-9903+0-0008. All 
specimens were prepared from spectrographically standardized materials. 


§1. INTRODUCTION 

N connection with an investigation of the susceptibility variation with 
if temperature of the transition metals it was necessary to determine the 

absolute susceptibilities of small specimens approximately 7-5 mm. long and 
4 mm. in diameter. With such small specimens it was advisable to adopt a 
comparative method rather than to attempt to determine the variation of field 
strength over a very limited range and thus evaluate the absolute force on a 
specimen in terms of susceptibility. 

A search of the literature showed that tantalum was a suitable metal for 
comparative work over the range required, especially as it is relatively cheap 
even in a highly purified state. This metal has an additional advantage as a 
standard in that it possesses a very small temperature coefficient of susceptibility. 

The only measurements on tantalum hitherto made appear to be those of 
Honda (1910) and Owen (1912) neither of whom was particularly concerned with 
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absolute accuracy. Their results are quite widely different, Honda giving 
y(18°c.)=0:93 x 10-§ and Owen y(18°c.)=0-81x10-®. It was, therefore, 
necessary to make an accurate investigation of the susceptibility of this metal 
at room temperature. 

§2. METHOD OF INVESTIGATION 

A rod of spectrographically standardized tantalum 10 cm. long and 4:5 mm. 
in diameter was obtained from Messrs. Johnson, Matthey Ltd. The major 
impurities are given in the Appendix. The dimensions of the rod made the 
Gouy method the obvious choice for an absolute determination. The field 
distribution of the available magnet had, therefore, to be investigated as accurately 
as possible. The method of field measurement adopted was a modification of 
that employed by Nettleton and Sugden (1939). Field strengths were obtained 
ultimately in terms of an inductance, an electric current, and the absolute 
dimensions of a search coil. All the electrical apparatus was standardized at the 
National Physical Laboratory at the conclusion of the work. ‘The measurements 
of search coil dimensions were made with standardized gauges. ‘The magnetic 
field was held constant by regulation of the magnet current which was maintained 
constant to one milliampere by means of an automatic regulator designed and 
built by Mr. J. C. Matthews. 

The remanent magnetism of the electromagnet used was rather high and an 
elaborate process of demagnetization was necessary between successive observa- 
tions to obtain reliable results. The experimental criterion of satisfactory 
demagnetization was that the residual field should be less than three gauss between 
the pole pieces. 

The field was measured along a symmetrical axis perpendicular to the line 
of the pole pieces for various excitation currents, the maximum being 5 amperes. 
Two search coils were used, it being found that the results thus obtained agreed 
to about one partinathousand. ‘The final results for the maximum field strengths 
were of the order of 10,000 gauss with a probable error (least squares) of + 3 gauss. 
The smallness of the probable error indicates the internal consistency of the 
measurements. 

§3. SUSCEPTIBILITY MEASUREMENTS 

The Gouy method employed for the susceptibility measurements requires 
no detailed description. Weighings were made in terms of calibrated weights 
standardized against N.P.L. certified weights. A Sartorius balance was employed, 
the ratio of the arms being determined. The tantalum rod was very slightly 
tapered towards the centre of its length, necessitating a small correction in 
evaluating the susceptibility from the weighings. No variation of susceptibility 
with field strength was observed, indicating absence of ferromagnetic impurities. 
The mean of twenty-two independent measurements gave the volume 
susceptibility as K x 10®= 14-171 + 0-009*. 

To convert this to mass susceptibility an accurate value of the density is 
required. ‘This was carefully measured using standard procedure and found to 
be p(20° c.) = 16-691 + 0-003 gm.cm-*, The value found by Keesom and Desirant 
(1941) is p(17°c.)=16-566+0-01 gm.cm-*. The wide divergence in these 
results is not understood: the higher value is in close agreement with that from 
x-ray data. 


* All errors quoted are by least squares, and where the result is dependent upon a previous 
result the error has been computed by a combination of errors in the standard manner. 
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With p=16-691 and atomic weight 180-88 the volume, mass and gramme 
atomic susceptibilities at 20°c., correcting by Honda’s (1910) result are: 
K x 10®°=14-171 + 0-009, x x 10°=0-8490 + 0-0006, x, x 10°=153-57+0-10. 


§4. THE MAGNETIC SUSCEPTIBILITIES OF PLATINUM, 
PALLADIUM AND RHODIUM 

Measurements were also made on small cylinders (0-75 cm. long, 0-4 cm. 
diameter) of platinum, palladium and rhodium which had previously been used 
in an extensive series of measurements on the low-temperature variation of 
susceptibility of these metals. Each of these had been prepared from spectro- 
graphically standardized material. The major impurities are given in an 
Appendix. Values of their susceptibilities at room temperature were obtained 
by a direct comparison with tantalum. Two cylinders were turned off the 
tantalum rod, one from the middle and one from the end, to the same dimensions 
as the cylinders of the metals to be compared. No difference in the magnetic 
behaviour of these two tantalum specimens was detected. 

Using a specially constructed holder the specimens were suspended in an 
inhomogeneous magnetic field and the force exerted on the specimen by the field 
was measured by direct weighing. A number of values of y at different 
temperatures was obtained and corrected to 20°c. using unpublished results 
on temperature variation previously obtained in these laboratories. The following 
results were obtained, all at 20°c.: platinum, x=(0-9712+0-0007) x 10-8, 
palladium, x =(5-231 + 0-004) x 10-6, rhodium, x =(0-9903 + 0-0008) x 10-§. For 
none of these metals was there any evidence of variation of susceptibility with 
field strength, and it may be concluded that there was no ferromagnetic con- 
tamination. ‘The palladium specimen was outgassed prior to the measurements 
by heating it zm vacuo at 100°c. for several hours and then at 450°c. for a few 
minutes, the glass enclosure then being sealed off and allowed to cool. The actual 
measurements were made as quickly as possible and there was no evidence of any 
absorption of hydrogen during the weighing. 


§5. COMPARISON WITH PREVIOUS WORK 

To facilitate the comparison with the values obtained by other observers 
a number of results are summarized in the Table. The figures in brackets are 
the values at 20°c. computed by taking the variation of susceptibility per 
centigrade degree as —0-0007, —0-010 and +0-0003 for platinum, palladium 
and rhodium respectively. 

It is seldom that any statement of accuracy is made and hence in the Table 
the significant figures quoted are the same as those given in the original papers. 
In view of the wide scatter of the results it would be unprofitable to attempt to 
assign limits of accuracy to the published values. 

Finke, Honda, Owen, Guthrie and Bourland all used torsional methods for 
susceptibility measurements, the first three making absolute measurements in 
terms of the field and its gradient. The determination of field gradient with 
a high degree of precision is difficult and consequently the order of agreement in 
the results of these investigators is not unsatisfactory. Guthrie and Bourland 
calibrated in terms of the susceptibility of water. 

Onnes and Oosterhuis investigated the variation of susceptibility of platinum 
and palladium with temperature using what is essentially a weighing method. 


They state that in the case of platinum it was difficult to determine the susceptibility 
Y-2 
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Susceptibility « 10° 
- Observer Pt Pd Rh 
Finke (1910) 1-054 5-478 1-031 
Honda (1910) 1-097 (18° c.) 5-79-(18° c.) 1-14 (18° c.) 
(1-096) (5:77) (1-14) 
Owen (1912) 0-80* (18° c.) 5-20 (18° c.) 1-08 (18° c.) 
(0°80) (5-18) (1-08) 
Onnes and Oosterhuis 
(1913 and 1914) 0-973: (17 -27c.) Sco dee) -= 
(0-971) (5:28) — 
Féex (1921) 1-035,.(17 1° ¢:) 5-264 (18-5° c.) — 
(1-031) (5-249) ae 
Guthrie and Bourland 0-982 (25° c.) Se 5x( 256063) 1-08* (25° c3) 
(1931) (0-986) (5-20) (1-08) 
Wucher (1950) — 5-338t (14:5° c.) — 
= (5-283) oN 
we 5-334 (13-6° c.) = 
we (5-270) ae 


Hoare and Walling (1950) 0°971"(20" 'c.) 5-23, (20° c.) 0-990 (20° c.) 


* Extrapolated to infinite field. + Specimen from Marret-Bonnin. 
accurately owing to its small value, and for palladium the curve of temperature 
variation showed more irregularities than would have been expected considering 
the accuracy of the experiments. 

Foex and Wucher each used a translational balance and standardized in terms 
of the susceptibility of solid manganese pyrophosphate which in turn was 
determined by comparison with solutions of manganese chloride of which the 
susceptibility was found by the ascension method. Whilst the two specimens 
of palladium used by Wucher gave values in satisfactory agreement at 20°Cc., 
the values differed by some 2% at 500°c, suggesting a physical or chemical 
difference in the specimens. 

It is remarkable that for the specimens used in the above investigations, 
excluding the present work, in all but four cases the source of supply was the firm 
of Heraeus. ‘The exceptions are one of Wucher’s specimens, supplied by Marret- 
Bonnin, and those of Guthrie and Bourland which were supplied by the Bureau 
of Standards. ‘The specimens used by Féex are stated to have been specially 
purified and free from iron. In the majority of the cases the susceptibility was 
investigated with various field strengths and the results of Owen (Pt) and Guthrie 
and Bourland (Rh) corrected to infinite field; no variation with field strength 
is noted in any other case where such an investigation is specifically mentioned. 
In no case is a spectrographic examination of the materials mentioned. Whilst 
it is not thought that the divergence in the results of later investigations is wholly 
attributable to impurities, it is considered that in work of this kind a spectrographic 
examination of the specimens used is well worth while. 
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APPENDIX 
Purity of Specimens used 


The following abstract of reports and private communications on the 
spectrographically standardized material used lists only those impurities whose 
spectrum lines are clearly visible. 

Tantalum. Columbium: 4100-923, 4079-729, 4058-938; nickel: all sensitive 
lines; tungsten: 4294-614, 4008-753. Columbium is the principal impurity. 
Tests on tantalum having the same source of origin might reasonably be quoted 
as indicating a typical and very approximate columbium content of 0:05%. - 

Platinum. No clearly visible lines. A quantitative assessment of some 
elements is given as follows: iron <0-001%, lead <0-001%, gold <0-001%, 
palladium <0-001%. 

Palladium. Silver: 3382-891, 3280-683, strong lines. Total impurities not 
more than 0-002%. 

Rhodium. Magnesium: 2852-129, 2802-695, 2795-53, copper: 3273-962, 
3247-540. Total impurities not more than 0-002%. 
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Rectifying Junction 
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ABSTRACT. If due allowance is made for the lowering by the image force of the height of 
the potential barrier in contact rectifiers, current-voltage characteristics can be calculated 
which are in good agreement with observed characteristics. On application of a voltage in 
the blocking direction, electrons (or holes) pass over the barrier owing to their thermal 
energy, the number of successful passes, i.e. the leakage current, being exponentially 
dependent on temperature. As the applied voltage is increased and the barrier lowered, the 
power loss corresponding to the increased leakage current will eventually raise the temper- 
ature in the thin barrier layer sufficiently to increase the current in turn until the cooling 
arrangement provided cannot cope with the increased power loss and the rectifier becomes 
thermally unstable. The conditions for this instability to arise are investigated quantita- 
tively. Amongst other parameters such as mode of operation, cooling, etc., the conductivity 
of the semiconductor, and in particular the height of the potential barrier, have a paramount 

_ effect on the maximum voltage that can be withstood. The simultaneous requirements for 
large resistivity and barrier height severely restrict the choice of materials which can be 
used for efficient rectifiers. 


N the present development of contact rectifiers considerable efforts are being 

spent on purifying the semiconductor material used. The physical properties 

of semiconductors are rather sensitive to traces of impurities, but particularly 
so are the properties of the barrier layer which fundamentally affect the character- 
istics of a rectifying junction. Germanium as used for transistors or for high-back- 
voltage rectifiers is a case in point. This need for high-purity materials has long 
been felt but has not so far been substantiated by any theoretical argument. 

According to the band theory of solids a potential barrier is set up at the 
junction of any n-type semiconductor* in contact with any metal of higher work 
function, and at the junction of any p-type semiconductor* in contact with any 
metal of lower work function (Angello 1949). In spite of the large number of 
metal-semiconductor combinations which should thus produce a barrier, practical 
contact rectifiers capable of rectifying reasonably large currents are confined to 
very few types; a somewhat larger number of metal-semiconductor combinations 
shows a certain amount of rectification when contact is made over a rather small 
area (point contact rectifiers). 

With the recent development of the theory of rectification it has become 
possible to discuss this problem quantitatively. It has been shown in the diffusion 
theory of rectification (Billig 1951) that the height of the potential barrier V, at 


* An n-type semiconductor has ‘ donor levels’ (easily ionized impurity centres) slightly below 
the conduction band, and conducts by free electrons, whilst a p-type semiconductor, with 
“acceptor levels’ above the full band, conducts by ‘ holes’. 
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the rectifying junction is reduced by the image force on the application of an 
inverse voltage — U to a value 


PEV TPIS ae aor ee aes (1) 
where Cocsixfel Vasa pee = a es (2) 


is the critical voltage, application of which would reduce the barrier completely. 
Here, e is the electronic charge, « the permittivity of the semiconductor and 
N the volume density of fixed space charges in the barrier region which is assumed 
to be uniform. Equation (2). holds for a space charge determined barrier 
(Schottky) to which the quantitative discussion below has been limited. Very 
similar results have been derived for a fully insulating (Mott) type of barrier. 

The power lost per unit area by the flow of leakage current on application 
of the inverse voltage is 


P=JU(U/V,)?* exp (— V/W) Lon: (3) 
BN ae . W me Var UN 
where J = . V4 (- Fe rm) a aeahe hi bree (4) 


and the ‘equivalent temperature potential’ W=kT/e. The power loss heats 
the rectifier and the leakage current rises correspondingly. This process continues 
and may lead to one of two results: (i) The cooling provided is sufficient to 
dissipate the heat generated and the rectifier reaches a certain steady temperature 
at which it may continue to operate indefinitely. (ii) The power lost at any 
temperature exceeds the heat that can be dissipated at that temperature; leakage 
current and temperature rise in turn until they become excessive and permanent 
damage is done to the semiconductor or the counter electrode. For a circuit 
element such as a contact rectifier, in which the power loss varies exponentially 
with the temperature, the criterion for thermal stability sets a limit to the maximum 
permissible temperature, the losses that can be dissipated at that temperature 
depending on the cooling arrangement. 

The highest inverse voltage — U, that can be sustained indefinitely is obviously 
of great importance and the effects of various parameters on its magnitude are 
discussed below. For a particular, fairly representative, type of selenium plate 
rectifier operating under A.c., a maximum voltage — U) =40 volts was calculated 
(Billig 1951) which is in good agreement with practical experience. ‘This value 
U, depends on various parameters, such as the operating conditions, type of load, 
ambient temperature, cooling arrangement, etc. In order not to confuse the 
issue we shall vary only one parameter at a time, for example when discussing 
the effect of the electrical conductivity of the semiconductor we maintain all 
other parameters entering the equations constant, in particular the height of the 
potential barrier, volume density of fixed space charges, heat dissipation, etc. 
This procedure is reasonable: potential barriers found at the junction of various 
combinations of semiconductors and metal electrodes are all of the same order 
of magnitude, say approximately 0-5 volt, and similarly space charge densities 
are of the order of 10! to 1018 cm’. 

We thus compare two rectifiers of identical size and construction, say of the 
plate type, but using different semiconductor materials: the electron mobility 
and hence the electrical conductivity o, of one, say selenium, is low, about 


2x 10-5 ohm-1 cm, and that, o, of the other, say germanium, rather high*, 
about 2 ohm-!cm-!. According to equations (3) and (4) the power lost in 
the germanium plate rectifier under the same voltage would be proportionately 
higher, by a factor 10°. If the selenium rectifier were stressed up to its limit 
(— U,=40volts,—j,=23 ma/cm?), the germanium rectifier under the same 
voltage, but with a leakage current 10° times higher, would immediately 
break down. To determine the maximum voltage U which the germanium 
rectifier could sustain indefinitely we maintain the critical temperature and thus 
the power loss fixed (at 0-92 watt/cm?). This leads to the transcendental equation 


(U/U,)"8 exp [a(U¥4— U,"4)]=0,/0, where «=V4/WUM4. ...... (5) 


This has been solved for the type of rectifier considered abovet, and the 
result is plotted in Figure 1 for a range of o from 2 to 2x10-* ohm™‘cm™. 
The voltage of the germanium plate rectifier would have to be reduced to 
1-73 volts for continued safe operation. 

Change in the type of metal used as counter electrode would affect the height 
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of the barrier Vy. Again maintaining the power loss as well as all other parameters 
constant, we have this time to solve 


(U/U,)"8 exp [8(U4— Up4)] =(Vqu/Va2exp (Va—Vay)|W] «+0-(6) 
where p= Val U."4 =(e/«)94(87N)"4, This has again been done for the selenium 
rectifier considered (o, =2 x 10-> ohm~ cm~) and the result is plotted in Figure 2 
showing how strongly the limiting voltage is affected by a change in the height 


* The large increase in conductivity is mainly due to the much higher mobility of the current 
er stipes rae than to any large increase in their volume density. 

e values of the various parameters are 0o,=2%10°ohm-!cm7}, V3 =0°527 volt 

Uz, =268 volts, U,=40-2 volts, a=5-21. : pe ; 
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of the barrier.. Actually the maximum temperature at which the rectifier can 
operate continuously is not independent of the barrier height, hence W in 
equation (6) is not quite constant; allowance for the corresponding change in 
the permissible power loss could easily be made, but the effect on the limiting 
voltage is relatively unimportant and has been neglected here. 

The way to overcome the difficulty of low stability inherent in the use of rather 
low resistivity semiconductors is to increase the effectiveness of the cooling 
arrangement. In practice this is done by reducing the contact area so that the 
power loss per unit contact area which can be dissipated by the rectifier is greatly 
increased. In the standard design of point contact rectifiers the cooling is in 
effect improved to such an extent that in spite of the greatly reduced contact area 
(~ 10° cm?) a leakage current of about the same magnitude as in a selenium 
rectifier of 1 cm? contact area can be tolerated, the area being reduced in abtou 
the same ratio as the conductivity of the germanium is increased. The maximum 
voltage that can be withstood before instability sets in is even higher in some 
germanium point rectifiers than in selenium plate rectifiers. Turn-over voltages 
exceeding 100 v. have been obtained. 
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ABSTRACT. It is shown that the self magnetic field of a high current discharge plays 
an important part in determining the radial electron density distribution. The effect 
of this magnetic field can be treated by adding a ‘ magnetic potential’ term Vn=fi wH dr 
to the electrostatic potential V, in the Boltzmann equation for the radial electron density 
distribution, w being the longitudinal electron drift velocity. The drift of electrons and 
ions to the tube walls is equalized at low arc currents by a radial electric field. Since the 
magnetic field reduces the drift of electrons to the wall whilst having little influence upon 
the ions, increasing the arc current leads to a progressive reduction of the radial electric 
field. Probe measurements in an arc at low pressures show that at currents over 
50 amperes, the radial electric field is zero except for a narrow sheath at the wall. The 
axial density is then mainly determined by the self magnetic field. A. theoretical 
expression is derived for the equilibrium density distribution assuming a balance between 
the diffusion rates due to the concentration gradient and the magnetic forces. Experiments 
in a mercury vapour arc at high currents support this theory. 


§1. INTRODUCTION 
YMBOLS USED: Suffix 1 refers to the electron, suffix 2 to the positive ion 
with m denoting mass, e charge, » number density, 7’ temperature, 
V temperature in equivalent units of potential; R is Boltzmann’s constant, 
#, longitudinal electron drift velocity, V, space potential, H, azimuthal magnetic 
field, V,, magnetic potential |’ #,H, dr, E radial electric field —(d V,/dr), F, radial 
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force on an electron, V,, potential difference across wall sheath, fo, pi, P2 neutral 
gas, electron and positive ion pressure, J, total arc current, J,,, positive ion 
current density at the walls, A total cross-sectional area in unit volume of normal 
atoms obstructing diffusion. Unless otherwise stated all equations are in 
electromagnetic units. ‘The temperature of the ions and the electrons will be 
defined by 4mc?=3kT/2, where c? is the mean square random velocity. This 
kinetic theory definition of temperature is applicable whether or not the gas is 
in a steady state (Chapman and Cowling 1939), but coincides with the 
thermodynamic definition when the gas is in equilibrium. 


Previous experimental investigations (Killian 1930) on the radial electron 
density and potential distribution in the positive column of an arc discharge have 
been confined to currents of less than 10 amperes. ‘Tonks (1939) has attempted 
to calculate the density distribution in high current arcs using methods of 
conventional diffusion theory. However, at low pressures (i.e. less than 
0-1 mm.Hg) where the positive ions fall radially towards the wall without making 
more than a few collisions with gas molecules, it is scarcely legitimate to assume 
that the ions have a Maxwellian velocity distribution. With increasing current, 
a diminution in the radial electric field may be expected, but until the energy 
gained per free path in this field is comparable with the thermal energies of the 
positive ions, diffusion methods are not applicable. 

Because the electron drift velocity to the walls is at least an order of magnitude 
less than the thermal velocity it is possible to assume that the approach to the 
equilibrium density distribution is very close, and therefore Boltzmann’s equation 
may be used to assess the importance of different potential fields, e.g. electrostatic, 
magnetic, gravitational, etc., in determining the electron density distribution. 

The general equilibrium radial charge density distribution follows from the 
equation g 


dps 
< (mRT;) + ( a Le (1) 


For example, at low pressures and currents, the only significant contribution to 
the (dp/dr) 5.14 term 1s that of the radial electric field. For T, constant, equation (1) 


then becomes RT, duldr—n.¢, fh =O, er (2) 
Integration leads immediately to the Boltzmann relation: 
ty =MoeXp(—eV,/RIy) | i eee (3) 


where V, is the electrostatic potential difference between two regions of density 
n, and M49. 

Experimental results of Killian (1930) with a current of 5 amperes confirm 
this dependence of the charge density on the radial variation of potential. The 
correction introduced by considering the magnetic field of the arc current is 
negligible for this value of current. An example of the low current distribution 
is illustrated in the experimental section. 

In this investigation it is assumed that the average potential energy of an 
electron increases towards the walls, due firstly to the electrostatic potential VE 
and secondly to a magnetic potential V,,. The magnetic potential arises fi6it 
the interaction of the longitudinal flow of the electrons and the magnetic field 
of the arc current. At high currents it is shown experimentally that the variation 
of V, across the greater part of the discharge cross section is negligible so that V,, 
is of major importance in controlling the density distribution. 
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§2. BOLTZMANN’S EQUATION WITH ELECTRIC 
AND MAGNETIC POTENTIALS 
The contribution of the magnetic field to the (dp/dr),,14 term can be determined 
- by considering the effective ‘ pressure’ it exerts upon the electrons. In cylindrical 
coordinates with the current flow parallel to the axis, the radial force F,, acting 


on a charge e,, is 
ee ee), dive, eae ree (4) 


where w,, is its longitudinal component of velocity, FE the radial electric field, 
H,, the azimuthal magnetic field. The pressure dp, on an element of volume of 
thickness dr distant r from the axis containing m, charges per unit volume is 


dp,=F, dr = +m Ee, —1,e,H,m, eeeree (5) 


where w, is now the average longitudinal drift velocity of the electrons at radius r. 
This is equivalent to defining w, by the equation J, =7,e,w, where J, is the drift 
current density in the discharge. Equation (1) now becomes: 

Ldn 147, | 4E  %H,w, 


ae Tar RT. kT, 
This mathematical procedure, though convenient, tends to obscure the 
fundamental processes. ‘The equilibrium state is attained when the net rate of 
flow across any surface due to diffusion is zero. Thus the terms in equation (6) 
may be interpreted as follows: the first two terms are proportional to the velocity 
of diffusion of charges due to a concentration gradient dn,/dr and a temperature 
gradient d7,/dr and the third and fourth terms to the rate of transport of charges 
in the electric and magnetic fields. Integration of (6) leads to 


=) ab Anes (6) 


n T re,Edr rt e, Hw, dr 
] ged Y ] a = | 1 = iste eak 
Ce oe aie.) eee 


Ny Tro 
where m,=m) when V,, H,=0, and T,=Ty. The ‘magnetic potential’ is 


defined by 


Vie i Hideaki oor (8) 
When 7, = 7, everywhere (7) becomes 
le ary Bs ails (9 
ae (V+ Va) RT. agen (a) 
Equation (9) can be written: 
ida sate eel O 
2 <exp| - gh-(V.+V a) |. 0) 


Equation (10) has been formulated without any reference to either the positive 
ions or the neutral gas atoms and it is necessary to enquire under what 
conditions it is valid. 

If the positive ions and the gas molecules are both in equilibrium under any 
potential fields which may be acting, the rate of flow of momentum across unit 
area is the same in opposite directions for each component. ‘Thus no net 
momentum can be communicated to the electrons and no additional pressure 
gradient can arise. However, in most practical discharge tubes, ions are flowing 
to the walls, and neutral gas returning from the walls into the discharge space. 
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Since a neutral gas atom usually makes many collisions with the wall before it 
collides with an electron, it is just as likely to deliver momentum to an electron in 
one direction as in any other. It follows therefore that, unless the mean free path 

of a gas atom for collision with an electron is of the order of the tube radius, this 
factor may be ignored. ‘The same argument applies to the positive ions drifting 
to the walls, viz. the free path for collision of a positive ion with an electron must 
be of the order of the tube radius before there is any tendency to deliver radial 
momentum to the electrons. At very high electron densities such as occur in 
the cathode spot of an arc or in a spark discharge, these factors may be important; 


§ 3. Feu LIER ELECTRON DENSITY DISTRIBUTION 
AT HIGH CURRENTS 

An equilibrium density distribution would result if there were no generation 
of charges and the boundaries were perfectly reflecting for both electrons and 
positive ions. If the net flux of particles is small compared with the random 
flux, the density distribution is, to a first approximation, an equilibrium one. 
This assumption is indeed made in the treatment of diffusion problems. A 
second approximation can be carried out to take into account a net flow of particles. 
In this section the equilibrium density distribution will be calculated assuming 
only a magnetic potential, as it is shown experimentally that the radial electric 
field is negligible across a large part of the discharge cross section. 

The equations to be solved are 


ty = Nig OkD.(—e,V {RI ,) ie) ye eee (11) 
vue | Birdie | (ASI yee ee (12) 
0 
2d r 
eo tee ae Newt GF. GM. eee (13) 


Measurements (to be published) on the axial electron drift velocity in a 
mercury arc as a function of the distance from the tube axis show that at pressures 
below 10-? mm.Hg and currents up to 150 amperes (the highest used) the drift 
velocity and electron temperature are independent of radius. This result is to 
be expected when the mean free path for electron scattering is greater than the 
tube radius. 

By differentiating (12) and (13), it is possible to eliminate H and to obtain 
a relation between nm andr. This is: 


d (r dn, 

a(x 2) + gmt =050)~ or) Peg aees (14) 
where g =(47e,?w,?/RT,). Writing x=Inr/r, and y =n,r? in (14) and integrating 
we obtain 1 dy\2 

€ Z) +2qy =constant. 


When r=0, y=0 and (1/y)(dy/dx)=2, so that the value of the constant is 4. 
Separation of variables and integration a gives 


tn eee lia | 
eee ere i MPa? ete, (15) 


where p is a constant. With m,=m) at r=0 this becomes ,=M49/(1 + pr?/r92)? 
where P= QMyo7 92/8 = 7e 12H ,2Nyo792/2RT}. Rasen LO} 
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The total arc current J, is 
To 
Jy =2n |" N,e,W4r dr =71€,W Nop ai(1 +p). 


Thus the axial electron density expressed in terms of the arc current is 
Ja 


Linea eric dn duvsfodien 4 RTh a abaaetnesaea (17) 
Eliminating 7,9 from p in equation (16) yields the final solution: 
ny 4d 61 » y2 7-2 
10 : ti (ee chk 4d 4 €W, 5 | Se OE ree (1 8) 


The average electron density 7, is simply m,)/(1 +A) where 
A — J 3 610,/(RT, — tJ 2610). 


In practical units, "eect (18) is: 


Jw 72-2 
1 1 
iB -[1+ Gas Toe El sariald 


where J, is in amperes, w, in cm/sec., V, in volts, and r in centimetres. 1/74 is 
plotted as a function of r for various values of the parameter A in Figure 1. 
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A comparison between the theoretical and experimentally determined electron 
density distribution is shown in Figure 2. The scales have been adjusted so that 
Equation (17) in practical units is 

Bee LO Og (20) 

M10 = F741 —5 x 10 ,0,)V 3) ot oi 
It is interesting to note that as J,w,—>2 x 10°V,, ny implying that the arc 
tends to contract into a narrow streamer at the tube axis. This behaviour may 
be compared with that termed by Tonks the ‘pinch effect’. With increasing 
arc current, 7, is generally observed to decrease whilst w, remains constant. 
If the discharge is to remain stable, this tendency must change in the neighbourhood 
of the critical current. If the degree of ionization approaches 100% the drift 
velocity is determined mainly by the positive ion scattering and both 7, and Wy 


(110) exp = (1210) theors 
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increase with increasing current (Mohler 1938) and it is not possible to decide 
whether the discharge will remain stable without a more detailed analysis. 


§4. THE RADIAL STEADY STATE MOMENTUM BALANCE 

Inastrictly equilibrium state the net rate of flow of momentum across any area 
is zero. In a plasma, however, ions and electrons are continuously generated 
and, in the absence of appreciable volume recombination, recombine on the tube 
walls. Consequently, in the steady state there is a radial flow of momentum. 
The ratio of the momentum carried by the positive ions to that carried by the 
electrons is m,/m,, since the drift velocities are equal (ambipolar drift). The: 
contribution of the electrons will therefore be neglected. 

The potential difference across the space charge sheath at the walls is. 
V,,=4V,1n(m.V,/m,V,). Thus even if V,=V,, the ions gain a greater part 
of their radial momentum in the wall sheath (V,,=6V, for mercury ions). 

At gas pressures where the free path for ion—gas atom collisions is small 
compared with the tube radius, the ions transfer their momentum to the neutral 
gas and set up a radial gas pressure gradient. ‘The momentum gained by the 
positive ions is then transferred to the walls by the neutral atoms. In this case 
it has been shown (Steenbeck 1939) that the difference in gas pressure between 
the walls and the tube axis is given by po(walls) — po(axis) =,)RT,. In the low 
pressure case, the rate of transfer of momentum by the positive ions is 7,m,0,”, 
and the positive ion wall current density is J,,,=m,€,v,. Thus the resultant 
pressure on the walls is p,(walls) =m.J,,02/e. The drift velocity of the ions at 
the walls is approximately v, =(2e,V,,/m,)"?. Thus the electron pressure at the 
plasma boundary can be expressed in terms of the electron temperature and the 
positive ion wall current: 


; . 
yg RT = matam (==) pas (- “ma ae (21) 


9 Mey 2 


Substituting for V,, from above and writing RT, =e,V, equation (21) becomes 


ho= VIB 6, (= In a) ur ui telvent Sepceenes . (22) 
Inserting the constants for mercury and putting V,=V, this becomes, with 
Z;, in milliamps, typ =3-2x 1007, (V2. = ee (22 a) 


At high currents, the electrons experience an inward force due to the self-magnetic 
field whereas the ions are little affected, and the axial electron pressure is greater 
than that given by equation (21). With n,=m,, at r=ry equation (18) gives 


110 =I, (mar) E + me 37 .Wy ik 
1 


eyVy — 3d W, 
I; 2 x 10°V. ; 
=a 11 _—2w 1 
x 10 Vp? E x oT | ’ eeeeee (23) 


for mercury ions, and with J,, in milliamps, J, in amperes, V, in volts, and 
w, in cm.sec"}. 
§5. APPARATUS 
Measurements were carried out in a Pyrex tube about one metre long and 
54 mm. internal diameter (Figure 3). A hollow graphite anode about 4 in. long 
and 1} in. diameter was used in most of the experiments, but it was later replaced 
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by a water-cooled nickel-plated copper anode A to enable the tube to be run at 
higher arc currents. <A’ is a small nickel anode used as ‘keep alive’, through 
which a constant current of about 2 amp. is passed. The cathode C consists 
of a mercury pool through the surface of which projects a #-in. diameter nickel 
cylinder used to fix the cathode spot. Probe readings are found to vary erratically 
when the spot is allowed to wander. A conical nickel shield B prevents direct 
streams of mercury vapour from the cathode spot entering the main length of 
the tube. The arc is started by applying a Tesla coil to a tungsten wire sealed 
through the tube T. 


Horizontal 
Movable 
Probe 


Inside Tube Diameter 
Z 54mm 


Ground Surface 


End of Tungsten Probe 
Closely Sleeved 


Probe to Mercury 
Pump 


Mercury Pool 
Cathode C {Water 


Figure 3. 


The spherical cathode compartment and half the neck are immersed in a 
vessel of water which is continuously stirred. The temperature of the water 
can be adjusted to any value in the range 0 to 40°c. by means of a copper cooling 
coil through which refrigerant liquid is passed from a thermostatically controlled 
refrigerator unit. 

The density of the neutral gas in the tube was calculated from the vapour 
pressure of the condensed mercury corresponding to the temperature of the 
water. Ifthe mean free path of the neutral gas is small compared with the diameter 
of the aperture between the cathode compartment and the tube, the pressure at 
the wall in the tube is equal to the vapour pressure of the condensed mercury, 
and the neutral gas density at the wall near the probe is 7/7 multiplied by the 
density corresponding to the condensed mercury temperature. If, however, this 
mean free path is comparable with or greater than the diameter of the aperture, 
thermal effusion will lead to a difference of pressure in the two regions and the 
pressure in the tube will be a factor (7y/7)'? higher than that in the bulb, and 
the density therefore a factor (7'3/Tyw)'* lower than that in the bulb. A heater 
wire is wound around the side arm P, so that when necessary its temperature can 
be raised above room temperature to prevent mercury from condensing in it. 
The electric field in the column is calculated from the difference in the floating 
potentials of the two probes P, and P, which are spaced a known distance 
(11-85mm.) apart. The plane collecting surfaces of the two probes are on the 
axis of the tube. A simple arrangement can be used for backing off the difference 
between their floating potentials against the voltage across a potentiometer, 
with a galvanometer to show equality. ‘The electric field can therefore be obtained 
at any moment from the reading of the voltmeter across the potentiometer. 
This was checked by taking the complete probe characteristics of P,; and P, 
and deriving the difference in space potentials at the two points in the usual 
manner. A movable probe P; is used to measure electron densities and space 
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potentials at any point in the cross section of the column. The collecting surface 
of P,is plane. The probe is made by sleeving a 14. mm. tungsten rod with vitreous 
silica tubing, the end of which was drawn down so that, with the rod in position, 
the end could be carefully ground plane with carborundum powder (Figure 3 (a)). 
In this manner, the collecting area of the probe is accurately defined except for a 
small error arising from the finite thickness of the positive ion sheath on the 
glass wall which tends to screen the probe from electrons approaching at small 
angles to the plane of the probe. Over the range of conditions used in these 
experiments, the positive ion wall currents were of the order 0-5 to 5-0 ma/cm?, 
and from Child’s space charge equation it is calculated that the sheath thickness 
ranged’ from about 0-17mm. to 004mm. The shielding will therefore be 
negligible. A probe P, in the same plane as the wall consists of a nickel disc 
0-95cm. in diameter and is used for measuring the positive ion current density 
at the wall. 

The initial outgassing of the apparatus consisted of running the discharge 
tube at a high current for several hours until the electrical characteristics were 
steady. The tube was subsequently kept continuously pumped, using a two-stage 
mercury pump. Inspection with a small spectroscope showed no impurity lines 
except sodium vapour at high currents. In the early experiments, the tube was 
pumped through a ?-in. diameter tube at the cathode end of the discharge tube 
and the rate of pumping limited the attainable arc current. Moreover, it was 
found that impurities tended to accumulate at the anode region due to the pumping 
action of the longitudinal electron drift current. The pumping tube was therefore 
removed to the anode end of the discharge tube (see Figure 3) and its diameter 
increased to 14 in. 

§6. MEASUREMENTS 

It will be seen from equations (12) and (13) that the magnetic potential V,, 
is directly proportional to the square of the axial drift velocity w, and it is therefore 
important to evaluate this parameter as accurately as possible, if a reasonable 
agreement with theory is to be expected. The most direct method would appear 
to be the use of the equation (11) for the total arc current J, in terms of the electron 
density distribution across the arc and the drift velocity w, although this is 
dependent upon the accuracy of the absolute value of peak electron density. 
Complete probe characteristics of voltage and current were taken at about 
ten positions traversing the arc. Correction for positive ion current was made 
when necessary and the plots on semi-logarithmic graph paper were linear over 
a wide range. Values of electron density , and space potential V, were obtained 
in the usual manner from these curves. The drift velocity could then be obtained 
from these data by graphical integration. These values of drift velocity were found 
to be in reasonable agreement with values obtained using the Langevin mobility 

Ss 
expression y= An. eX itp 
(2m,kT,)1 A 
where A is the averaged cross section for electron collision assuming a Maxwellian 
velocity distribution corresponding to an electron temperature 7, and X is the 
longitudinal electric field. This averaged cross section was computed using 
data from Brode (1929) for the variation of the cross section of electrons of uniform 


energy and by integrating over the Maxwellian velocity distribution as described 
by Killian (1930). 


Self Magnetic Field in High Current Gas Discharges 353 


The results of the probe characteristics taken at arc currents of 5 amp. 
and 50amp. are shown in Figures 4 and 5 respectively. It will be seen from 
Figure 4 that, the space potential V,, plotted in units of V, where V,=RT;/e,, 
is of the expected parabolic shape, whereas in Figure 5, for the high arc current, 
the space potential has become uniform over the arc except at the walls. The 
lack of symmetry in the curves about the axis is probably due to the self-shielding 
of the probe as it approaches the far wall of the tube. Equation (9) may be 
written In(m,/7,9)=|(V,+V,)/V,|. When the magnetic term is zero, the space 
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potential curve should be coincident with a curve of In(m/79). ‘This is shown 
with good agreement in Figure 4 for an arc current of 5 amperes where the 
magnetic potential term is very small (two points at the walls with the computed 
magnetic correction are shown). In Figure 5 (50 amperes arc current) there is a 
wide deviation in the In(m,/n,)) plots, but a plot of (V,+V)/V, brings the 
agreement considerably closer. If a slightly higher drift velocity had been 
assumed, an even better agreement with the theory would have been attained. 
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Figure 6. Comparison of measured axial electron density with that 
computed from equations (20), (22 a) and (23). 


In Figure 6 the variation of axial electron density as measured by a probe 
is shown as a function of arc current. It departs considerably from the linear 
relation predicted by low current theory (‘Tonks and Langmuir 1929). Much 
better agreement is obtained on the basis of the magnetic potential theory using 


equation (20), and this is further improved by using equation (23) which 
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incorporates a value of electron density at the wall deduced from radial momentum 
considerations. Figure 6 also includes a plot of the electron density at the plasma 
boundary (equation (22 a)) for reference. 


§7. DISCUSSION 

As the are current increases, the magnetic field becomes more and more 
important in reducing the radial drift of electrons and the radial electric field 
which maintains ambipolar drift at low currents decreases until it is effectively 
zero. Concentration gradient diffusion is then responsible for ambipolar drift 
to the walls. A detectable reversal of the radial electric field (and thus a negative 
arc core) cannot in fact occur whilst the ions have temperatures comparable with 
that of the tube walls: a reversal of 0-04 volt would be sufficient to bring about 
a considerable reduction in the radial drift velocity of ions having an equivalent 
temperature of 100° c. 

The absence of the radial electric field leads to an increase in the lifetime of a 
positive ion with the result that electron-ion collisions and collisions of the 
second kind may increase its energy. If this is so, the ion temperature. may be 
expected to approach the electron temperature at sufficiently high arc currents. 
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The Relationship between the Radial Concentration Variations 
of excited Atoms and Electrons in a Discharge in Thermal 
Equilibrium 
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ABSTRACT. 'The excitation and Saha equations applicable to gaseous electrical 
discharges in thermal equilibrium are used to obtain convenient expressions for calculating 
the radial concentration variations of electrons and excited atoms when the radial intensity 


variation of one spectral line is known, and when the temperature on the axis of the discharge 
can be estimated. 
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§1. INTRODUCTION 
N most cases electron concentrations in gaseous electrical discharges are 
| obtained as mean values of the whole of the discharge cross section. For a 
more detailed study of such problems as Stark effect, the effect of positive ions 
on electron mobility, and current densities, or self absorption, a knowledge of the 
spatial variation of electron concentration is required (see for example Schulz 
1947, Edels to be published). The measurement of such a spatial variation is 
by no means simple, although it may be obtained for example from Stark effect 
investigations on radially symmetric discharges if the radial intensity variation of 
the spectral line under investigation is known (Grimley and Saxe, to be published, . 
Hopwood, Craig and Craggs, to be published). However, if such an intensity 
variation has been obtained it is a relatively simple matter to derive either the radial 
intensity variation of another spectral line, or the radial variation of electron 
concentration if the discharge is assumed to exist in thermal equilibrium. 

The assumption of thermal equilibrium in a gaseous electrical discharge in 
effect postulates that the energy gained by the charge carriers from the electrostatic 
field between the electrodes is passed on by collision to the neutral particles in the 
discharge to such an extent that the energy of all types of particle is distributed 
according to the Maxwell—Boltzmann law with the same temperature, that is with 
the same mean kinetic energy. Since energy transport is continually taking place 
in the discharge, thermal equilibrium cannot be complete, and in effect a difference 
between the temperature of the electrons and ions and the other particles must 
exist. However, for many discharges with a high gas density there is considerable 
theoretical and experimental evidence to show that this temperature difference is 
small, and in such cases thermal equilibrium may be assumed as a close approxi- 
mation. With such a discharge the radial energy transport can also be assumed to 
have only a small effect on the energy distribution of the particles, so that in each 
small volume of the discharge the Maxwell—Boltzmann distribution is maintained 
at a temperature characteristic of the spatial position of the considered volume. 
The radial energy transport is then associated with a radial temperature distribution 
to which the formulae of thermal equilibrium may be applied. For a further 
discussion of this problem and a description of the theory and experiments 
associated with the subject reference may be made to a report on the temperatures. 
of an electrical discharge in a gas and to the references cited there (Edels 1950). 
In a discharge in which thermal equilibrium is assumed, both the radial spectral 
intensity and electron concentration variations are functions of the common 
radial temperature distribution, through which they may be related. The 
relationships so obtained are developed in the following section. 


§2. RADIAL RELATIONSHIPS IN A DISCHARGE IN 
THERMAL EQUILIBRIUM 
We consider an atom having energy levels s, t, i (the ionization level), such that 
the magnitudes of the energies to the ground state p are given by hy,, << hy, < hy. 
We assume now that the intensity of radiation at a radius 7 obtained by transition 
from level s to some lower state is known as a fraction of its intensity at the discharge 
axis and is given by F. We wish then to obtain the radial variation of radiation 
intensity f of a transition from level t, and the electron concentration variation g as 
functions of F. It will be appreciated that the radial variation of the level 


populations will also be given by F and f. 
Zao 
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For a discharge in thermal equilibrium the intensity of radiation of the frequency 

Vp at r is given by the excitation equation as 
- | —hy 
In= Nowe Av gAgy XP f ai | sks Oe (1) 

where N, is the concentration of atoms in the ground state, g, and g, are the 
statistical weights of the states s and p, A,, is Einstein’s probability of spontaneous 
transition between the two states, and 7, is the temperature at 7. Ornstein 
and Brinkman (1934) derive and discuss the validity of this equation and its 
applications to arc discharges. 

If the temperature at the discharge axis is taken as 7, from equation (1), then 
the radial variations of intensity are given by 


ai To —hrg (To . ] 1 
[i= Th {exp RT, (7 1 eed ee coe’ oe ( a) 


pss 7! {exp Sr" (7-1) |}. rea (18) 


Now Saha has shown (1920, 1921, 1923) that for a gaseous volume at 7 in thermal 
equilibrium, a relationship between the electron concentration V,, concentrations 
of positive ions NV, and neutral atoms NV, and the temperature may be obtained. 
This relationship, developed by Fowler (1929), is 


NN, ge(2am)3? oes ae » | Fr hy; | 
a a Ae osha 2 
ei OS re (2) 
where g, and m are the statistical weight and mass of the snes and U,(T) and 

U,(T) are functions containing the statistical weights and partition terms of the 
energy states of the ions and atoms respectively. From this equation, asso 
a neutral plasma (IV, = NV ,), we can derive the function g as 


g= (Zz) {exp ES (7 - 1) |. ETS (2 a) 


The derivation and validity of equations (1) and (2) is given by Edels (1950). 

The derivations of equations (1a), (1b) and (2a) assume a constant pressure 
discharge which is applicable in the majority of cases. For the early period of a 
spark discharge or for other special sources in which constant density conditions 
may exist, the factor outside the exponential is different. However, as will be seen 
later, this term is not of primary importance and the main results will remain valid. 

Thermal equilibrium, as stated previously, can normally only be obtained in 
high density discharges which have high temperatures associated with them. 
Under these conditions it is found experimentally (Brinkman 1937, Kruithof 1943) 
that the temperature falls a relatively small amount over the whole of the visible 
radius. ‘This is indeed to be expected if hv/kT,~10, since the exponential terms in 
(1a), (16) and (2) control the functions. With this condition it is possible as a 
first approximation to take the functions F, f and g as solely proportional to the 
exponential terms. This gives 


fot = ( 2) ee ee ee (3) 


where «=14,/v) or v;/2v,, for f and g respectively. In such gases as hydrogen, 
helium, etc., the electron concentration variation will differ markedly from the 
intensity variation, since for these cases v;/2v,,~4, so that if the visible limit is 
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taken as F~0-01 theng~0-1. This effect may appreciably affect measurements 
of average current densities if the temperature does not fall off rapidly beyond 
the visible radius. 

Assuming always the validity of equations (1) and (2), the possibility of estimating 
the error in (3) and making more accurate calculations depends upon the value of 
hy./R To, and for a given spectral line upon the accuracy of the knowledge of Ty. 
It is evident that if 7 is known, equation (1) determines, at every value of F, the 
value of 7)/T,, and hence of f and g. We can improve upon the accuracy of the 
relationship (3) by making further approximations or by utilizing the exact 
relationships between the radial variations. These relationships are given by 


fe Ote sa (T Lyre he eras (4) 


where the values of « are as above and B=v,/Vsp and (4+ »4/2v.p) for f and g 
respectively. 

Now the variations of F with 7,/T) for different values of hv,,/RT) can be 
determined from equation (1a) and the results plotted in the form of curves as shown 
inthe Figure. Thus if a function F has been obtained for a line of some value of 


: N Tr/To=1 a\2 
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hy,,/RT), the curves can be used to obtain the corresponding values of T/T) and 
hence from equation (4) the functions f and g may be obtained. 

Now with T)~5,000 °x., RT)~0-5 ev., so that for the high ionization potential 
elements hy,,/kTy>~20. An examination of the Figure shows that for this case the 
intensity has dropped to 1% when 7,,/T, has reached only 0-8. Since also B = 
V4g[Vsp~0-2 it is clear that for these elements equation (3) for fis applicable. For the 
electron concentration case B~0-9 and an error is introduced if equation (3) is used. 
The error is not great since it is a percentage error of g and reduces as g increases. 
Thus with g~50%, the value given by equation (3) will be ~52%. For the low 
ionization potential elements however hyv,,/kT, may approach the value5. Inthis 
case the curves show T,,/Ty~0-6 when F=5%. In this case B =1,/v,, may still be 
small, in which case the approximation for f, given by equation (3), may be used. 
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For the function g however the true relationship given by (4) must be used. In 
* general it may be stated that with values of hy,,/RT) greater than or equal to 10 
the simple relationship given in equation (3) may be used without undue error. 

More accurate approximations to the true relationship can be obtained by 
determining 7,,/T, in terms of F from equation (1 a) by the method of iteration and 
substituting the eSNg approximations in panes (4). 

Putting 7,/T,.=-7 and the first approximation 7 , = 1 we obtain for subsequent 
approximations 


F,=(1—ln-F)=s 
TF, =x+EIMT ny 
where £=RT,/hv,,. This gives the general mth approximation as 
fy OF By =(F{x+Eln7F,, 1)? and f/f, or g/g2={T,/T.(1 —£1n F)¥. 
From equation (1 a) we find also that 


(1—éIn F)=(1)/T,)(1 pg In(T)/T,,)}, 


ee ae TL, To\" 
so that Sor] =41- ee F zh bare (6) 


Now for this case for F=0-05 when hy,,/RT)=5 ; f/f, or g/g. ={0-94}", and it is 
clear that only a small error will result if the second approximation 1s used, and this 


ae fa OF fa =(F)*/{1—(RT [Ryn FY. see 7) 


This approximation may therefore be used for all cases when hyv,,/RT, >5, and will 
thus cover most elements. 

To summarize, an accurate plot of f or g can be obtained from equation (4) in 
conjunction with the curves given. For high ionization potential elements with 
hy,,/RT,) greater than or equal to 10 the first approximation f or g=(F)* 
may be used. For values 10>h»,,/kRT,)>5 the second approximation 
for g=(F)*/{1 —(RT,/hy,,) In F}’ should be used, although the first approximation 
may still be used for the f values. 

The analysis shows that in all but extreme cases the electron concentration does 
not fall off radially as rapidly as the radiation, so that care must be taken if intensity 
variations are being used as an indicator of electron concentration variations. 
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The Interpretation of Radiosonde Data in Relation to 
Cosmic-Ray Intensity Variation 


Nicolson and Sarabhai (1948) have recently discussed difficulties in interpreting 
quantitatively the semi-diurnal variation of cosmic-ray intensity in terms of the theories of 
atmospheric oscillations. Duperier (1948) has also discussed the possibility of relating 
vertical movements of the higher levels of the atmosphere to the ‘temperature effect’ 
in cosmic-ray intensity in particular. : 

One major difficulty in this latter study lies in correct interpretation of the radiosonde 
results. It is difficult to disentangle true heating and resulting expansion of the free air 
from any positive temperature error due to insolation of the thermometer element itself 
when the ascents take place in daylight. 

The ascents from two radiosonde stations, Downham Market (524° N., $° E.) and 
Larkhill (513° N., 2° W.) are here analysed for the year 1948. The ascents take place at 
03, 09, 15 and 21 hours G.m.T. daily. The ascents at 03 and 21 hours occur in darkness 
and the 09 and 15 hours in sunlight throughout the year. The two latter are liable to 
thermometer error. 

The data are taken from the Daily Weather Reports (Upper Air Section) issued by the 
Meteorological Office. The heights at fixed pressure levels are computed from temperatures 
and pressures measured by the Kew radiosonde (Dymond 1947). 

In the Table are given : Column (1). H, the seasonal mean height of the pressure levels 
300 mb. and 100 mb. at 03 and 21 hours G.m.T. (The nocturnal change in height from 
21 to 03 hours is small compared with the apparent daylight changes, and it is convenient 
to consider the former ascents together. The height given in the first column will 
approximate to conditions at midnight.) Column (2). AH,, the apparent increase in height 
of these two pressure levels since midnight at 09 hours G.m.T. Column (3). AH», the 
apparent increase in height at 15 hours G.mM.T. Column (4). f,, the percentage of the total 
apparent height change which occurs above 300 mb. at 09 hours. Column (5). fe, the 
corresponding value at 15 hours. 


(1) (2) (3) (4) (5) (6) 
Season H(km.) AH, (m.) AH, (m.) f:(%) f5(%) serch nbs) 
Spring 9-11; +3 +13 90 78 300 
(Feb., Mar., Apr.) 16-16, +29 +58 100 
Summer 0-235 =, +20 aS 77 300 
(May, June, July) 16°17, +68 +85 100 
Autumn 9-31, +14 +23 78 69 300 
(Aug., Sep., Oct.) 16°38, +64 +73 100 
Winter 9-075 +9 +9 72 56 300 
(Nov., Dec., Jan.) 16-03, +32 +27 100 
Standard deviation of AH values: c. +60 metres at 300 mb. ; +120 metres at 100 mb. 
Probable error of seasonal AH values: c. +4 metres at 300 mb. ; --8 metres at 100 mb. 


It will be noticed that the major part, 75° on the average, of the apparent daylight 
increase in height of the 100 mb. pressure level takes place above the 300 mb. pressure level. 
It is just in these higher levels that there is least certainty about the magnitude of whatever 
temperature error there may be due to insolation on the thermometer assembly of the 
radiosonde. If in daylight the thermometer were to read too high between 300 mb. and 
100 mb., by about 3° x. in winter, rising to about 1°. in summer, then the values AH, 
and AH, for 100 mb. would be reduced to the values given for 300 mb. _ Insolation errors 
approaching this order of magnitude are not impossible. Furthermore, it would be 
extremely difficult to account for the measured daylight temperature increases. ‘The flux 
of solar energy in the lower stratosphere is insufficient to provide the apparent diurnal 
heating by absorption of radiation, unless our knowledge of the composition of the 
atmosphere at these levels is very false, and the dynamic processes required to produce 
the observed diurnal changes would be extremely complex, when it is remembered that 
diurnal surface pressure changes must remain within the small limits observed in our 
latitudes. There is at the moment little direct evidence for the existence of such diurnal 
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dynamic effects as are required, and also little experimental assessment of the magnitude of 
radiation errors in actual flights. (The meteorological problem is discussed in more detail 
in a paper to be published shortly in the Quarterly Journal of the Royal Meteorological 
Society.) 

For ascents which take place symmetrically about local noon, and which have probably 
similar radiation errors, it appears that the differential heating, which ought to be real, is 
a slowly increasing function of height. ((AH,-AH,) has an average for the whole year 
in this analysis of 4 metres at 300 mb. and 10 metres at 100 mb.) If these figures should 
correspond more closely to the real variation in AH with height than do the measured AH 
values, then correlations of AH with simultaneous ground temperature changes will 
necessitate the choice of higher levels of the atmosphere to attain the coefficient of 18 metres 
per centigrade degree arrived at by Duperier (1948) from discussion of diurnal changes. 

It should be noticed also that, whereas the apparent diurnal temperature changes in 
the lower stratosphere are rather closely symmetrical about local noon, the maximum of 
the ground level diurnal temperature change is shifted to about 15 hours at all seasons. 


Conclusions. 


(1) In the absence of evidence concerning the reality of the apparent diurnal changes 
in height of the isobaric levels in the lower stratosphere, it is best to disregard these changes 
in correlating H with ground temperature to arrive at a value for the ‘temperature effect’ 
in cosmic-ray intensity. 

(2) Correlation of annual variation in the monthly mean value of the height of the 
100 mb. pressure level with ground temperature would seem to be fully justified. The 
radiosonde ascents which have taken place since 1938 bear out as substantially correct 
Blackett’s assumption (1938) of a seasonal change in the monthly mean value of the height 
of this pressure level of the order of 500 metres, accompanying the mean summer—winter 
variation in ground level temperature of about 10° c. 

(3) Since day-to-day changes of the height of the 100 mb. pressure level can be of 
the order of 300 metres (i.e. approaching the magnitude of the seasonal change in monthly 
mean height), it would seem worth while to correlate these directly with changes in cosmic- 
ray intensity. 

(4) The meteorological data are now so abundant—four ascents daily from six stations 
in the British Isles—that there is no difficulty in synchronizing cosmic-ray observations 
with meteorological ascents. Since ground level temperature changes are so closely 
dependent upon local conditions, it is perhaps better to refer the ‘temperature effect’ to 
larger fractions of the atmosphere (see Millican and Loughridge 1948). Data are given by 
the Meteorological Office at intervals of 50 mb. from ground level to 200 mb. (approxi- 
mately 14 km.) at 20 mb. intervals to 100 mb. (16 km.), and thereafter at 10 mb. intervals 
to the top of the ascent. 

(5) The measured apparent diurnal changes in the height of the 100 mb. pressure level 
do not explain either the magnitude or phase of the diurnal variation in cosmic-ray intensity. 


The Clarendon Laboratory, Rabe Kaye 
Oxford. 
3rd February 1951. 
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Properties of Amorphous Selenium in the Infra-Red 


The optical constants of amorphous selenium in the visible region are well known 
(Barnard 1930, Becker and Schaper 1944), but only qualitative results in the infra-red 
seem to have been recorded (Becker and Schaper 1944, Elliott, Ambrose and Temple 1948). 
The purpose of this note is to give some results of recent measurements. Specimens of 
amorphous selenium in the form of prisms and parallel-sided slabs were made by allowing 
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pure melted selenium to solidify rapidly between two optically worked glass plates. The 
Figure shows the transmission T of a typical parallel-sided specimen as a function of wave- 
length. On the long-wave side of 0-7 transmission is high and is determined in this 
thickness (0-162 cm.) almost entirely by the reflection at the two surfaces. 
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The absorption coefficient « at any one wavelength may be determined from the slope 
of the (log 7, thickness) curve, the error due to multiple reflections being small. The 
lowest values of ~ found in measurements of a number of prisms are : 


A (u) 0-8 1-0 2-0 4-0 
a (cm-) 3-3 1:9 0-55 0-52 


Some specimens gave higher values, probably owing to partial crystallization. 

Refractive indices for prisms were obtained by measuring the angle of minimum 
deviation at 0-819 and 1-014, using an infra-red image convertor tube. The mean values 
obtained from prisms of angles between 5° and 15° are for A=0°819p,, n=2:589+0-005, 
and for A=1-014y, n=2-520+0-005. ‘The index values were extended to longer wave- 
lengths by an interference method using thin films. At wavelengths beyond 2 the index 
attains a nearly constant value of about 2°45. The low value of absorption suggests the 
use of amorphous selenium as an optical material for infra-red wavelengths. It has the 
advantage of ease of preparation and is stable in the atmosphere. 

Thanks are due to Prof. R. W. Ditchburn for helpful advice and to Messrs. Standard 
Telecommunication Laboratories, Ltd., for a grant (E. W.S.). One of us is a member of 
the Royal Naval Scientific Service (H. A. G.). We also wish to thank Prof. Lark Horovitz 
and Dr. K. Frank, with whom we had discussions during the 1950 Reading Conference 
on the Properties of Semiconducting Materials. 
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A Note on the Dielectric Dispersion in Polycrystalline 
Materials 


Dielectric dispersion has recently been observed in certain polycrystalline solids, 
particularly semiconductors, due to their inhomogeneous structure. Thus, Chasmar 
(1948) has observed a large decrease in the resistance of PbS and PbSe films with increasing 
frequency, and interpreted this as a bridging of the high resistance between the intergranular 
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boundaries by their capacitive reactance. Similarly, Volger (1951) has explained the 
decrease of both capacitance and resistance of certain ceramic materials by their grain — 
structure. German work along similar lines on selenium has been reported (F.1.A.T. 
Report 1946). 

To explain such results, a model is assumed (Figure 1 (a)) consisting of low-resistance 
grains G (resistance Rg) in series with boundaries B of much higher resistance Rp shunted 
by acapacitance Cy. The high value of Rp may be due to the presence of an air gap reducing 
the area of contact between grains as postulated by Volger, or of a layer of material of higher 
resistivity, for instance amorphous selenium in the case of microcrystalline selenium. 


1g 3 
R Cc 
T, C, 
(c) 
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Figure 1. 


At high frequencies the boundary resistance is shunted by its capacitance so that only the 
grain resistance is measured, whilst at low frequencies resistance and capacitance of the 
boundaries are measured. It is clear that an abnormally high permittivity is obtained 
if the large capacitance value measured at low frequencies is referred to the total thickness 
of the sample. 

The model can be elaborated (Figure 1 (b)) by including the capacitance Cy of the 
grains, which becomes important at higher frequencies (Cg< Cp, owing to the much smaller 
thickness ot the boundaries). It can then be shown that the apparent permittivity and 
loss angle vary with the frequency in a similar way as in the case of Debye dispersion. - 

The expressions resulting for the equivalent capacitance C and shunt resistance R of the 
model circuit shown in Figure 1 (6) are 


1+ w?R,?Cp? 1+w?R2C,C 
R=h) —— 2 —— C=Cp—— 5 
1+ eR RbCp 1+ oRECy 
AE 1 1+ w*RgRpCp? 


RCo RyCyw 1+e°RZCCh 


under the simplifying assumptions Rp > Rg and Cp > Cg. 
Owing to the power loss component (wCR))~1 due to the D.c. conductivity, the value of 
tan 6 first drops as the frequency is increased, goes through a minimum 


tan 8min=2(Rg/Rp)? at @min=1/Cp(RpRg)”, 
and then passes through a maximum value 

tan 8max=3(Cp/Cg)'? at ®max = 1/Rg(CgCp)*?. 
The ‘ a.c. component’ of tan 4, i.e. 

inoievbosst oRgCp 
RCw | RypCw  1+0*R2C,Cy’ 
is of the standard Debye form. 
Figure 1 (e), (d), (f) illustrate the variation with frequency of apparent resistivity p’, 


permittivity «’ and power factor tan8 for a typical polycrystalline semiconductor with 
parameters assumed as given in the Table.* 


* These values correspond in order of magnitude to standard selenium rectifiers. 
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Polycrystalline Se per cm? area 


Crystalline Amorphous 
grains boundaries 
Permittivity K 10 10 
Resistivity p (ohm cm.) 10? 10? 
Total thickness t (cm.) 10-2 Os 
Resistance R (ohms) 1 < 10? 
Capacitance C (pF.) 10? <— 108 
Time constant 7 (sec.) 1o=9 <n. ne 
High Low 
frequencies frequencies 
Apparent 
permittivity x’ 10 108 
Apparent 
resistivity p’ (ohm cm.) 10? 10? 


Another case of abnormally high dielectric constant has recently been reported (Laurila 
1950) for thin ‘sputtered’ layers of CaF,. The apparent dielectric constant when measured 
shortly after deposition of the layers was very high (550 to 4,400) but, within a day or so, 
dropped to more normal values (5 to 10). At the same time, however, the a.c. and D.c. 
resistance of the samples increased by a factor varying from 15 to over 700 for different samples. 
Such behaviour might be explained by assuming single crystallites of high resistivity 
embedded in a matrix of comparatively low resistivity which practically short-circuits the 
gaps between individual crystallites. Such a low-resistivity matrix may simply consist of 
condensed moisture or of some other conducting impurity deposited between the grains. 
As this is driven off or is replaced by continued growth of the. crystallites the overall 
resistance would be expected to rise and the capacitance to drop to more normal values. 

It should be mentioned that the model of an inhomogeneous semiconductor differs 
only in some detail from the model of an inhomogeneous dielectric (the Maxwell—Wagner 
mechanism of dielectric loss), which also gives rise to a set of Debye equations. 

Thus from dielectric measurements alone it cannot be decided what mechanism is 
causing dielectric dispersion, but often the magnitudes involved will point to the above 
model as the most likely explanation of a particularly high dielectric constant. ‘The true 
electrical parameters of the grain material are then measured only at high frequencies; 
but useful information on the relative thickness and resistance of the boundaries may also 
be deduced from the extreme values of the tan8 curve over the significant range of frequencies. 
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A Surface Effect in the Creep Behaviour of 
Polycrystalline Lead 


There is much experimental evidence to show that the mechanical properties of metal 
single crystals are influenced by the condition of the surface, in particular by its freedom 
or otherwise from oxide layers (Andrade 1949, Andrade and Randall 1948). In poly- 
crystalline wires a surface effect of a somewhat different kind is to be anticipated. It is 
tolerably clear now that the f, or transient, flow and the «, or permanent, flow are, in the 
main, due to two different physical processes, although they are not independent phenomena 
(Andrade 1948), the f flow being particularly associated with glide within the crystal 
grains and the « flow with intergranular movement. Now glide within the grains is, of 
course, influenced by surface conditions, and in a polycrystalline material the grains at 
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the free metal surface are in a privileged position, in that one boundary is free, while those 
in the neighbourhood of the surface are somewhat less constrained than those well in the _ 
interior of the metal. Hence a surface effect on creep is to be anticipated in the sense that 
the surface layers should offer relatively small resistance to creep. Such an effect has been 
established. 

The effect will clearly be easier to detect if only transient flow is present. The creep 
of commercial lead at room temperature is closely represented by the equation 

l= BE expi(ct)) a) a eee ees (1) 

with finite values of x, and so is the creep of the purest lead, but here recrystallization under 
stress takes place (Andrade 1948 b). It has been found that the addition of 0-05% 
tellurium to pure lead gives a metal which, at room temperature, shows no permanent flow 
except at the highest stresses, the creep being, in general, accurately represented by (1) 
with x=0. The present experiments have been carried out with such lead. Certain 
pecularities have been observed at very low stresses, and at the highest stresses « has a 
finite value. The following results refer to a range of stresses (about 50 to 95 kg/cm?) 
within which f varies with stress in the s gmoid way described in the early papers (Andrade 
1910, 1914) and since confirmed by Sinclair and by Jolliffe working in the Carey Foster 
Laboratory (to be published shortly). 

Experiments were carried out with wires of four different aber ters, viz. 1, 2, 3 and 
5 mm. approximately, the grain diameter being very nearly the scme, 0-2 mm., in each case. 
Figure 1 (a) shows the values of 8 at different stresses for the four wires at a temperature 
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Figure 1. (a) f against stress for different diameters; (b) 8 against effective stress, according to 
formula (2), for different diameters. 


of 31:5°c.: the discussion of the peculiar behaviour at low stresses we reserve for later 
full publication.. The stresses are, of course, the average longitudinal stresses for the 
whole cross section. It will be seen that to produce the same transient creep in different 
wires requires an average stress which is higher the thicker the wire, which is what we 
should expect if there were a relatively weak boundary layer. 

The results can be represented by assuming that the average stress S required to give a 
chosen value of f, i.e. a chosen transient creep, is given by 


S=S974/7=70) 0 eee eee (2) 
where a is the average grain diameter, 7 the radius of the wire, S, the value of S for very 
large values of r, when the surface effect is negligible, and m is a constant to be found, na 
representing in an arbitrary way the equivalent thickness of the weak boundary layer. 
Taking B=2:50x10-*sec-"8, for example, for which, from low to high diameter, the 
corresponding measured stresses are 66:0, 72:25, 75:20 and 77-3 kg/cm?2, the value n=0-301 
gives for Sp in the four cases 82:13, 82-19, 81-81 and 81-74 kg/cm?, which vary from the 
average by less than 0-3%, a constancy as satisfactory as can be expected. 

Figure 1(b) shows the values of f for all the observations, for the four different diameters, 
plotted against the effective stress as deduced from (2) with n=0-301. It will be seen that, 
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neglecting the abnormality represented by the broken lines (which in the range of effective 
stress 47 to 110 kg/cm® applies to the 1 mm. diameter wire only and will, as has been said, 
be discussed in full publication) all the observations for effective stresses above 50 kg/cm? 
are very satisfactorily represented. The value of 7 is also a reasonable one : the surface 
effect can be represented by supposing that a layer of thickness about one-third of a grain 
diameter offers no resistance to creep. This is clearly an arbitrary way of representing 
the surface effect: m may be called ‘the effective depth coefficient’ of the surface layer. 
It remains to be tested whether n is constant with different grain size. 

The non-uniformity of behaviour over the cross section has been established by x-ray 
examination, the method being to reduce the diameter of the strained wire by successive 
cold etching with 1 part hydrogen peroxide to 3 parts acetic acid (Worner and Worner’s 
solution) and take back-reflection photographs at each stage with characteristic copper 
radiation. For a specimen of initial diameter 1:95 mm., subjected to slow creep at 65° c., 
in which a total strain of 6:5°% was produced in 250 hours, two close well-resolved rings 
due to the copper Ka doublet are obtained at all depths, there being, however, a slight 
broadening of the lines in the case of the original surface. For a 3 mm. aiam. specimen more 
rapidly extended at room temperature (final strain 10°% attained in 10 minutes) the surface 
effect is the same, but the broadening of the lines increases progressively with increasing 
depth, the Ka doublet being resolvable at a depth of 0:16 mm. but not at a depth of 0:25 mm., 
the broadening being due, no doubt, to heavy internal stresses. At a depth of 0-78 mm. 
(final diameter 1 mm.) the smeared ring is very broad. With more rapid extension to 10% 
final strain the surface gives the same clearly resolvable rings as in the case of slow extension, 
but the interior exhibits recrystallization, evidenced by well-marked discrete spots in the 
tings. ‘This recrystallization is accompanied by internal stresses, evidenced by the general 
broadening of the rings. At higher rates of extension, however (10% final strain attained 
in 15 seconds), the interior shows large-scale recrystallization and little evidence of internal 
stresses, as would be expected, the stresses being relieved by the adjustments involved in 
recrystallization. This is illustrated in the series of photographs shown in Figure 2. 


(a) (b) (c) (d) (e) 
Figure 2. Back-reflection photographs of tellurium lead. 
(a) Surface of wire, unetched, after 10% extension in 15 sec.: diameter 2°73 mm. 

(6) Specimen etched down to 1°75 mm. diameter. (c) Etched down to 1°40 mm. diameter. 

(d) Etched down to 1:0 mm. diameter. (e) Untreated surface of another part of the original 

specimen after the previous photographs had been taken, as contro]. Reflection from (220) 

and (311) planes. 

That the difference between the internal photographs and the surface photographs is not 
due to general recrystallization taking place subsequent to the straining, and progressing 
during the time occupied in the preparation of the specimens by etching, is shown by the 
fact that surface photographs taken subsequent to the internal photographs with another— 
untreated—part of the specimen showed the same sharp rings as initially. 

The slow extension at 65° means a small extension due to f flow, the rapid extension 
at room temperature a practically complete f flow (Kennedy 1949). 

The x-ray work, which has established that a fairly sharp change in the type of structure 
takes place at a depth somewhere about 0-1 mm. below the surface, is continuing, with the 
object of establishing in more detail the nature of the transition. It is cited here to emphasize 
the real existence of a surface effect in the flow of polycrystalline metal. 
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Note added in proof. In a paper which has just reached us, T’. Nishishara and S. Taira 
(Mem. Fac. Engng. Kyoto Univ., 1950, 12, 90-118) report observations on the non-uniformity 
of stress across the cross section of steel and aluminium specimens deformed plastically in 
tension. As far as they go, their qualitative deductions are in agreement with what we here 
record. 
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ABSTRACTS FOR SECTION A 


The Shp, Twinning, Cohesion, Growth and Boundaries of Crystals, by H. WiLMAN. 


ABSTRACT. The nature and stability of intercrystalline boundaries are discussed and 
conclusions as to preferred relative orientations are illustrated by ball-bearing and bubble- 
raft models. 

Electron-diffraction evidence from crystals of layer-lattice, ionic, metallic and van der 
Waals (organic molecular) types is presented, showing that crystal pairs and symmetrical 
or unsymmetrical triplet or multiplet crystal groupings occur, having a common lattice 
row, but in relative azimuthal orientations at intervals agreeing with the predictions. ‘The 
component crystals were thus in coherent metastable contact of the type discussed, and in 
some cases strong secondary elastic scattering showed that the crystals were extensively 
superposed on the plane which was normal to the common row. 

The seven types of electron-diffraction pattern observed from such groupings, and 
the nature of the specimen preparation, suggest that the origin of these groupings is a 
mechanical deformation of the crystals during or after growth, by a process of ‘rotational 
slip’ which has not hitherto been explored or defined in nature. Such slip is directly 
demonstrated macroscopically as a deformation process in potassium ferrocyanide trihydrate 
and gypsum, the slip being produced on the planes parallel to the highly perfect cleavage. 

The relation between rotational slip and the known types of translational slip and 
dislocation theory is indicated, and in particular the ‘ deformation bands’ investigated 
especially by Barrett and his collaborators are concluded to be essentially rotational slip 
bands initially at low deformation, though modified by translational types of slip at higher 
deformation. 

A valuable new insight is also provided into many observations on crystal growth and 
properties, especially the nature, determining conditions, stability and disorientation of 
crystals growing epitaxially on single crystal substrates, and into the nature of the deforma- 
tion caused by unidirectional abrasion of crystals. 


The Variation with Temperature of the Electrical Properties of a Degenerate 
Electronic Semiconductor as exemplified by Cadmium Oxide, by R. W. WRIGHT. 


ABSTRACT. Strong strip specimens are prepared from CdO powder such that they 
possess an electrical conductivity of the order of 300 ohm cm™! at room temperature. 
Simultaneous measurements of the Hall constant R and the electrical conductivity o are 
made throughout the temperature range 0—500° c., and the thermoelectric power dE/dT 
is measured immediately afterwards on the same specimen. ‘The properties are found to 
resemble those of a metal : R is found to have a constant value independent of the tempera- 
ture. dE/dT is approximately 60 uv/deg. c. at room temperatures, whilst o decreases 
with temperature, its value at 500° c. being roughly half that at 0° c. 

The theory of degenerate semiconductors is developed, and since the values of R 
indicate that there is a constant number of free electrons in each specimen, this result is 
incorporated into the comparison of the theory with the experimental results. Satisfactory 
agreement is obtained between theory and experiment. Agreement is also obtained when 
the theory is compared with results obtained by another investigator for the thermoelectric 
power of CdO. The approximate form of the theory predicts that the thermoelectric power 
and the resistivity should be directly proportional to the absolute temperature. ‘The 
experimental results are found to confirm this. CdO behaves as an excess degenerate 
semiconductor, with properties strongly resembling those of a metal. 


Electron-Bombardment Conductivity of Dielectric Films, by F. ANSBACHER and 


W. EHRENBERG. 

ABSTRACT. Films of some dielectrics become conductive under the influence of 
electron bombardment. This conductivity is investigated both for steady and for 
alternating bombardment in its dependence particularly on temperature, current density 
and the potential difference across the dielectric. The gain, i.e. the ratio of current 
through the film to the bombarding current increases strongly with temperature. 
On arsenic sulphide films, the gain reaches 40,000. The effect is discussed in the light of 
theories of conduction in dielectrics and of recent work on allied subjects. 
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The Specific Heat of Liquid Helium at Temperatures between 0-6° and 1-6°K., 
by R. A. Hutt, K. R. WILKINSON AnD J. WILKS. 
ABSTRACT. A detailed account is given of experiments made to measure the specific 


heat of liquid helium between 0:6 and 1-6°k. Below 1-4° the specific heat may be expressed 
as C=0-024T °2 cal/gm/deg. Some implications of this result are briefly discussed. 


The Total Cross Section of Beryllium, Aluminium, Sulphur and Lead for Neutrons 
of Energies from 2 MeV. to 6 MeV., by G. H. STAFFORD. 

ABSTRACT. The total neutron cross sections of beryllium, aluminium, sulphur and 

lead have been measured in the range of energies from 2 Mev. to 4 Mey. and from 5 Mev. to 

6 Mev. by a transmission experiment. Resonances were observed in all four elements. 

The results obtained for lead and for beryllium are compared with those expected according 

to the one-body theory for light nuclei. 


New Radioactive Isotopes produced by Nuclear Photo-Disintegration, by 
F. D. S. BUTEMENT. 

ABSTRACT. Photonuclear (y, n) and (y, p) reactions produced by 23 Mev. x-rays from 

a synchrotron have been used to prepare new radioactive isotopes, the method being 

particularly advantageous for neutron-excess isotopes in the region Z=63 to 79. Whenever 

possible the isotopes were chemically identified. The following table summarizes the 


results. Blement 
bron Chemically Half-Life Assignment 
Identified 
Se As 9 minutes ?9As 
Sion Sm. 8 minutes 14391 OF ot 
Eu Eu 15 hours BI Dye 
Gd — 20 minutes 159Ru or Gd™ 
Dy — 14 minutes TOE Opec aD 
Er = 96 minutes Mens (HO OF -oo HO 
Yb —- 19 minutes eH gp cetteal Ip aaigopc mth-d bao) 
Hf Lu 22 minutes TS ore se 
W Ta 48 minutes hod Des} 
W ARES 116 hours teen 
Os: ea Re 17 minutes 189Re or 191Re 
Pt Ir 66 minutes Poo TTOne wht 
Hg Au 27 minutes 201Au or 2% Au 


Definite mass assignments, hitherto uncertain, have been made for 18-hour ™®Gd_ and 
48-minute ?9°Au. 


On the Proton Component of the Vertical Cosmic-Ray Beam at Sea Level, by 
M. G. Mytroi and J. G. WILson. 


ABSTRACT. Approximately 1% of the vertical cosmic-ray intensity at sea level is shown 
to consist of protons. These form a momentum spectrum of index —2°8 in the range 
10° mev/c< p<5 x 10% Mev/c. and have an absorption range in lead, 160+40 gm.cm-. 
Using existing measurements at high altitude, the attenuation in the atmosphere of protons 
in the spectral region investigated is found to take place with a range 140+-10 gm.cm-?. 

There is no negatively charged component of strong nuclear interaction in the vertical 
beam of intensity as great as 0:1°% of the whole beam. 


Further Measurements of the Charge Ratio of jx- Mesons at Sea Level, by B.G. OWEN 
and J. G. WILSON 


ABSTRACT. Extended measurements of the charge ratio (positive/negative) of near- 
vertical -mesons at sea level show that it increases with meson momentum from about 
1-17 at 10? Mev/c. up to about 1:26 at 4103 Mev/c. At higher momenta the statistical 
accuracy of the measured charge ratio is limited by the instrumental rate of collection, and 
it cannot yet be definitely decided whether the charge ratio decreases with increasing 
momentum, although the evidence now available favours a slow decrease. 

The general implications of the form of variation of charge ratio with momentum are 
indicated. 
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(a) | () 
Figure 1. 


(a) Wax impression of crater in lead block due to jet from tin cap. 
(6) Wax impression of crater in lead block due to jet from steel cone. 
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(4) 1,000 c/s. pulse- 


Figure 5. 20 msec. vibrator voltage pulse. 


a) 300 c/s., rod, mountings. (6) 600 c’s., rod mountings. 


C/S., Spring Mountings (d) 600 c’s., spring mountings 


1ons recorded at 130 ft., with the vibrator in various conditions. 
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Figure 6, Vibrations recorded at various frequencies and distances. 
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